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RESUMO 


EXISTÊNCIA E CONFIGURAÇÃO DE CICLOS LIMITES EM SISTEMAS 
DIFERENCIAIS PLANARES 


Neste trabalho, apresentamos resultados sobre existência de ciclos limites 
em sistemas diferenciais planares, configuração de tais ciclos e a dinâmica de 
sistemas não-lineares nas vizinhanças de órbitas periódicas. 
PALAVRAS-CHAVE: Existência de ciclos limites. Configuração de ciclos limites 
em campos vetoriais polinomiais planares. Teoria Qualitativa. 


ABSTRACT 


EXISTENCE AND CONFIGURATION OF LIMIT CYCLES IN PLANAR 
SYSTEMS 


In this dissertation, we present results on the existence of limit cycles in planar 
differential systems, configuration of such cycles and the dynamics of nonlinear 
systems near periodic orbits. 

KEYWORDS: Existence of limit cycles. Configuration of limit cycles in planar 
polynomial vector fields. Qualitative theory. 
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SUMÁRIO 


INTRODUÇÃO 


A teoria qualitativa das Equações Diferenciais Ordinárias (EDO's) busca entender 
o comportamento de soluções de EDO's sem explicitá-las. Um episódio muito importante 
no estudo dessa teoria se deu em 1900, em Paris, no Congresso Internacional de Ma- 
temática, quando o alemão David Hilbert enumerou, pela primeira vez, sua famosa lista 
de 23 problemas matemáticos, nos quais os mais proeminentes pesquisadores da época 
deveriam se ocupar para impulsionar o desenvolvimento da matemática como ciência. Até 
o momento, o XVI problema ainda permanece em aberto. 

A segunda parte do XVI Problema de Hilbert tem o seguinte enunciado: Qual o 
número máximo de ciclos limites de um sistema diferencial planar polinomial, e o que 
podemos dizer sobre as posições relativas dos seus ciclos limites? 

Embora esta teoria tenha avançado muito ao longo das últimas décadas, ainda 
existem questões de difícil compreensão sobre o assunto. Por exemplo, J. ÉCALLE [05] 
e Y. IPYASHENKO [16] mostraram no início da década de 1990, que o número de ciclos 
limites de um sistema polinomial planar é sempre finito. Resta ainda determinar uma cota 
superior para o número de ciclos limites de um campo polinomial planar que dependa 
somente do grau do polinômio que define o campo. 

O conceito de ciclo limite para campos vetoriais planares, introduzido por Poincaré, 
modela oscilações não lineares em processos físicos. No final da década de 1920, Van der 
Pol, Liénard e Andronov, estudando oscilações não lineares de fenômenos elétricos, anali- 
saram certas equações diferenciais ordinárias de segunda ordem verificando a ocorrência 
de ciclos limites. Após tal verificação, matemáticos e físicos estudaram extensivamente a 
existência e unicidade, entre outras propriedades, destes ciclos limites. 

O principal objetivo desta dissertação é demonstrar uma espécie de recíproca do 
XVI problema de Hilbert, provada por J. LLIBRE e G. RODRÍGUEZ [07] em 2004. Isto é, 
dado um conjunto finito qualquer de curvas de Jordan no plano, duas a duas disjuntas, é 
possível explicitar um sistema de EDO's polinomial que tem exatamente tais curvas como 
ciclos limites. Mais do que isso, o resultado estabelece uma relação entre o número de 
curvas de Jordan dadas e o grau do polinômio que define o sistema de EDO's. 

O método 'Averaging' também é abordado no presente texto. O 'Averaging' é uma 
ferramenta clássica que nos permite estudar a dinâmica dos sistemas diferenciais não- 
lineares sob forçamento periódico. Ele tem uma longa história começando com as obras 
clássicas de Lagrange e Laplace, que forneceu uma justificativa intuitiva do método. A 
primeira formalização desta teoria foi feita em 1928 por P. Fatou. As contribuições para a 
teoria do 'Averaging' foram feitas na década de 1930 por Bogoliubov-Krylov. 

O método do 'Averaging' é uma ferramenta para o estudo de sistemas não lineares 
com coeficientes periódicos que consiste essencialmente em determinar a existência e 


estabilidade de ciclos limites em um sistema de EDO's, a partir do estudo das singularidades 


Introdução 


de outro sistema, chamado de sistema promediado. Mais precisamente, a existência e 
estabilidade de cada singularidade do sistema promediado determina, sob certas hipóte- 
ses, a existência e estabilidade de órbitas periódicas do sistema original. Em particular, o 
método permite estudar a quantidade de ciclos limites de um sistema diferencial. 
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ASPECTOS GERAIS 


Neste capítulo estabelecemos a notação e resultados básicos que serão utilizados 
no decorrer do texto. 

Sejam O um subconjunto aberto de R x R" e A um subconjunto aberto de R”. Um 
ponto de R x R” será denotado por (t x), onde te Rcom x=(x, x, :::,x)E Re lx 
denota uma norma de Rº. 

Sejam k m, n=1,2,..., O espaço vetorial das aplicações X: A > R” que têm todas 
as derivadas parciais de ordem menor ou igual a k contínuas será denotado por CA, Rm), 
e definimos C”(A, R”) como o espaço vetorial de todas as aplicações X tais que XE CHA, 
R”), vk E (1, 2, ...>. Denotaremos por C“(A, R”) o espaço vetorial de todas as aplicações 
analíticas X: À > R”. Os elementos de CHA, R”) são chamados de aplicações de classe 
Cr. 

Será denotada por Df (x) a matriz jacobiana de f avaliada no ponto x € A. 


2.1 ORDEM DE CRESCIMENTO DE FUNÇÕES 
Utilizaremos as notações de O e o, como segue abaixo, para estimar certos termos 
em uma determinada equação. 
Definição 2.1.1 Sejam A c Rº aberto contendo o ponto0 ef g: A > Rº com g(x) * 
0, VxE A. Escrevemos 
(9) = O(g(x)) 
se existirem constantes Ce ô >0 tais que 


O<lxii<xôó<llf(MIl<cC, 
II g(x) | 


quando 


Observação 2.1.2 Vejamos algumas propriedades que decorrem diretamente da 
Definição 2.1.1: 


1. Se f(x) = A(g(x)) e (x) = O(g(x)), então f(x) + Ex) = O(g(9); 

2. Se f(x) = O(g (2) e Lx) = O(g(x)), então f09).(x) = O(g,(X).9,09); 
3. Se f(x) = O(lxl”) e0 <n=<m, então f (x) Axl” = O(llxil”r-"); 

4. Se f(x) = O(llxll”), então lixil"f (x) = O(llxil”+"); 

5. Se f(x) = O(llxll”), então f (x) = O(lixll"), para n< m; 
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Analogamente, valem propriedades semelhantes para o. Ainda, vale que se f(x) = 
o(Iixil”), então f(x) = O(llxll”). 


2.2 TEOREMAS DE JORDAN E GREEN 


O Teorema da Curva de Jordan desempenha papel importante ao longo deste texto 
e sua demonstração pode ser encontrada em J. MUNKRES [08]. 

Uma consequência útil é o Teorema de Schonflies, que afirma que a região limitada 
determinada por uma curva de Jordan é homeomorfa ao disco unitário no plano. A versão 
aqui apresentada pode ser encontrada em P. V. ARAÚJO [13]. 

Um instrumento fundamental de que necessitaremos é o Teorema de Green, tam- 
bém conhecido como Teorema de Stokes no plano. A versão apresentada aqui do Teorema 
de Green pode ser encontrada em M. G. SOARES [11]. 


2.2.1 Os teoremas da Curva de Jordan e Schonflies 


O Teorema da Curva de Jordan possui muitas versões; utilizaremos aqui a mais 
conveniente para os resultados que seguem. 

Teorema 2.2.1 (Teorema da Curva de Jordan) Seja y uma curva fechada simples em 
R2. Então, Rº — y tem exatamente duas componentes conexas, cuja fronteira comum é y. 

Como R? não é um conjunto limitado, uma das componentes de R? — y é limitada e a 
outra não. Chamaremos de Inty (interior de y) a componente limitada de R? — y, e de Exty 
(exterior de y) a componente ilimitada de Rº — y. 

Uma curva satisfazendo as hipóteses do Teorema 2.2.1 é usualmente chamada de 
curva de Jordan. 

Teorema 2.2.2 (Teorema de Schonflies) A região limitada por uma curva fechada 
simples é homeomorfa a um disco aberto. 

A demonstração pode ser consultada em P. V. ARAÚJO [13]. No Apêndice desse 
texto apresentamos uma versão do Teorema de Schonflies para grafos, pois é mais 
conveniente para os objetivos desse trabalho. 


2.2.2 Teorema de Green 


Para enunciarmos o Teorema de Green necessitamos do conceito de compa- 
tibilidade de orientação entre um subconjunto do plano e sua fronteira. 

Definição 2.2.3 Um subconjunto não vazio U c Rº chamado um domínio se U é 
aberto e se, dados dois pontos quaisquer p e q em U, existe um caminho de classe C'!, 
inteiramente contido em U, cujos pontos inicial e terminal são, respectivamente, p e q. 

Definição 2.2.4 Suponha que U c Rº é um domínio de um caminho de classe C', 
e seja Vc U um subconjunto fechado e limitado, cuja fronteira OV consiste de um número 
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finito de curvas de Jordan de classe C' e tal que VioV é um domínio. Para cada uma 
dessas curvas, adotamos o sentido de percurso para o qual o interior de V está sempre à 
esquerda quando a percorrermos. Nessas condições, dizemos que Ve dV têm orientação 
compatível. 

Definição 2.2.5 Se y : [a, b] > U é um caminho de classe C' em U, com U c Rº 
aberto, sendo y(? = (x(), Y(D) para te [a, b], definimos a integral de f ao longo de y (ou 
integral de linha de f ao longo de y) por: 


b 
[1 = far + vdy := ] [u(a(t), y(t)a'(t) + v(x(t), y(t))y'(t)]dt. 


Sey,*Y,*... *y, é a justaposição dos caminhos de classe C' em U, a integral de f 
ao longo desse caminho é definida por: 


[DR RAS RARA RE: 


Exemplo 2.2.6 Para o conjunto V da Figura 1.1, U=R2, dV=y,Uy,U y,. Para que 
a orientação seja compatível, y, deve ser percorrido no sentido anti-horário e y,, Yy, ambos 
no sentido horário. 

Teorema 2.2.7 (Teorema de Green) Sejam U c Rº um aberto conexo ef: U > Rº 
uma aplicação suave. Seja Vc U um subconjunto satisfazendo: (i) Vé fechado e limitado, 
(ii) a fronteira OV de V consiste de um número finito de curvas de Jordan de classe C!, OV 
=y,U...Uy, e (ii) VidV é um domínio. Suponha que Ve dV têm orientação compatível. 
Então, escrevendo f(x, y) = (u(x, y), V(x, Y)), temos que 


! = / ud + vdy = ) (5 Es 5) dady. 
av av |: dr Oy 


Figura 2.1 — Orientação das curvas do Exemplo 2.2.6. 


Fonte: Elaborado pelo autor. 
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2.3 CAMPO DE VETORES E FLUXO 


Detalhes dos resultados apresentados nesta e nas demais seções deste capítulo 
podem ser encontrados em J. SOTOMAYOR [09]. 

Definição 2.3.1 Um campo de vetores de classe C%, 1 <k=<o ouk=aw,é uma 
aplicação X: A > Rº de classe C*, onde A c R" é aberto. 

Ao campo de vetores X está associada a equação diferencial 


x = X(09 (2.1) 


e vice-versa. Aqui o ponto denota a derivada com respeito ao tempo t A partir de 
agora, nesta seção, X sempre denota um campo de vetores definido em um aberto A de R". 

Se x, € A é um ponto tal que X(x,) = 0, então x, é dito ponto singular, caso contrário 
é chamado ponto regular. 

As aplicações q E C'(I, A), onde [é um intervalo da reta, tais que 


(1) = X(o(t) 


para todo t E |, são soluções, trajetórias ou curvas integrais de X. A imagem de uma 
trajetória é denominada órbita. Denotaremos por y, a órbita passando pelo ponto p € A. 

Chama-se solução máxima de (2.1) a toda solução q definida num intervalo /, deno- 
minado intervalo máximo de q, tal que se y é outra solução no intervalo Jcom y :J> O 
talquelcJe q = WII, então |= Je, consequentemente, q = UV. 

Apresentamos a seguir o teorema de existência e diferenciabilidade de soluções 
máximas com relação às condições iniciais. 


Teorema 2.3.2 Seja X um campo de vetores de classe C&, k=1,2,...,0 ouk=0, 
definido num aberto A de Rº. 


1 (Existência e Unicidade de soluções máximas) Para cada x e À existe um 
intervalo aberto |, onde está definida a única solução máxima q ,: 1, > A, que 
satisfaz dy/dt= X(y), V(O) = x. 


IH (Propriedade de grupo) Sey= q(t sete l, então | =| -t=tr-trel) 
e qp(s y=d(t+s,x) para todo s E , 


Hl. (Diferenciabilidade em relação às condições iniciais) O conjunto D = «(t, 
x); xe A, te |) é aberto em R"*! e a aplicação q : D-> R" dada por q (t, x) é 
de classe C*. Mais ainda, q satisfaz a equação 


DD, v(t x) = DX(p(t x) º D,p(t x), 
para todo (t, x) E D. 
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Aqui o símbolo D, q (t x) representa a matriz n x n cujas colunas coincidem com as 
n últimas colunas de Dq (t, x). Já o símbolo DD, q (t, x) denota a derivação das entradas 
de D,q(t, x) com relação à variável t. A aplicação q : D => A dada no Teorema 2.3.2 é 
chamada de fluxo de X. 

Se qe Vi então existe t, E A tal que q= p(t, p), logo p(t q) = qp(t+t, pe : =, 
= |. Pelo Teorema 2.3.2 segue que y, = y,. Em outras palavras, dadas duas órbitas de X, 
então elas coincidem ou são disjuntas. 

Observação 2.3.3 Como qp(t, q(-t, x)) = x, segue que para cada t E |, fixo, a 
aplicação x > q(t x) é um difeomortismo. 

Seja q uma curva integral de X. Dizemos que ela é periódica se existe um número 
real T>0Otal que p(t+ T) = qy(?, para todo te Re q é não constante. Nesse caso, di- 
zemos que a órbita associada é uma órbita periódica. Observamos que uma singularidade 
não corresponde à uma orbita periódica. 

O próximo resultado assegura-nos que, dada uma solução qy de um campo de ve- 
tores X, ela é ou um ponto, ou uma linha homeomorfa à R ou uma órbita periódica. 


Teorema 2.3.4 Se y é uma curva integral máxima de (2.1) em |, verifica-se uma 
única das seguintes alternativas: 
q é uma bijeção; 


HH I=Regy éconstante; 


HM I=Req éperiódica, isto é, existe T>0 talque p(t+ T)= q(), vte R. 


O próximo resultado estabelece que órbitas com intervalo máximo limitado 'esca- 
pam” de compactos. Logo, toda órbita limitada está associada a uma solução cujo intervalo 


maximal é igual a R. 


Teorema 2.3.5 Sejam A c R” um aberto contendo um ponto x, X: A > Rº de classe 
C'e(a, 8) o intervalo máximo da solução x para o problema de valor inicial 


x = X(%), 
x(0) =x 


o 


Assuma que B <º. Então dado um compacto Kc A, existe te (a, B) tal que x(t) não 
pertence a K. 


Teorema 2.3.6 Se X: A > Rº é um campo de classe C!, definido em um aberto sim- 


plesmente conexo A c Rº, então as seguintes afirmações são equivalentes: 


1. existeH:A> R de classe C' tal que VH = X; 
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2. J,X=0 para todo caminho fechado simples em U; 


8. Ro ) X não depende do caminho que liga os pontos (x,, Y,) e (x, ). 


z0s40 


2.4 EQUIVALÊNCIA E CONJUGAÇÃO DE CAMPOS DE VETORES 


A fim de comparar retratos de fases de dois campos vetoriais, faremos algumas 
considerações. 

Definição 2.4.1 Sejam X, e X, campos definidos nos abertos U, e U, de Rn, 
respectivamente. Dizemos que X, e X, são topologicamente equivalentes se existe um 
homeomorfismo h : U, > U, que leva a órbita de X, em órbitas de X,, preservando a 
orientação. O homeomorfismo h é dito uma equivalência topológica entre X, e X,. Se h for 
um difeomorfismo de classe C', com r= 1,2, ...,º our= w, dizemos que X, e X, são C”- 
equivalentes. 

Mais precisamente, se pe U, e Yo é a órbita orientada de X, que passa por p, temos 


h(y,) = Up onde a, é a órbita orientada de X, que passa por h(p). 


Definição 242 Sejam X, e X, campos definidos nos abertos U, e U, de Rº, 
respectivamente. Sejam ainda q: D > Re q,:D, > R" os fluxos associados aos 
campos de vetores X, e X,, respectivamente. Dizemos que X, e X, são topologicamente 
conjugados se existe um homeomorfismo h: U, > U, tal que h(p (t x)) = q (t, h(x)), para 
todo (t, x) E D,. O homeomorfismo h é dito uma conjugação topológica entre X, e X,. Se h 
for um difeomorfismo de classe C”, dizemos que X, e X, são C'-conjugados e, neste casso, 
h é dito uma C"-conjugação entre X, e X,. 

Cada uma das duas definições anteriores estabelece uma relação de equivalência 
no conjunto dos campos de vetores definidos em U, e U,. Além disso, no caso da Definição 
2.4.2 temos necessariamente que trajetórias correspondentes têm o mesmo intervalo 
máximo. Uma equivalência h leva ponto singular em ponto singular e órbita periódica em 
órbita periódica. Além disso, se h for uma conjugação, o período das órbitas periódicas é 
preservado. Observe ainda que toda conjugação é uma equivalência, mas o contrário não 
é válido. 

Exemplo 2.4.3 Sejam X,=(y -x) eX,=(2y, 2x). Temos X, e X, são C'-equivalentes, 
pois basta tomar h : RP > Rº a função identidade dada por h(x) = x. Mas X, e X, não são 
conjugados, pois o período das órbitas difeomorfas não são iguais. 

A seguir apresentamos uma caracterização de conjugação diferenciável. 

Proposição 2.4.4 Sejam X, e X, campos vetoriais de classe C' definidos nos abertos 
U eU, de Rº, respectivamente, eh: U, > U, um difeomorfismo de classe C;, r=1,2,...,0 
ou w. Então h é uma conjugação entre X, e X, se, e somente se, 


para todo xe U,. 


Aspectos gerais 


2.5 FLUXO TUBULAR 


Nosso próximo objetivo é mostrar que todo ponto regular de um campo de vetores 
possui uma vizinhança na qual o campo é diferenciavelmente conjugado ao campo 
constante. Este resultado é bem conhecido e trata-se do Teorema do Fluxo Tubular. 

Definição 2.5.1 Sejam X: A > R” um campo de vetores de classe Cie Ac Rr! um 
aberto. Uma aplicação diferenciável f: A> A de classe C* chama-se seção transversal local 
de X (de classe C*) quando, para todo a E A, Df(a) e X(f(a)) geram o espaço R". Seja > = 
F(A). Se f: A> > é um homeomorfismo então dizemos que > é uma seção transversal de 
X. Quando 0 € A e escrevemos p = f (0) dizemos ainda que > é uma seção transversal de 
X emp. 


Observação 2.5.2 Sejam p E A um ponio regular de Xetv, v,*:*,V,, X(p)) uma 


n=1º 


base de R". Seja g: Rr! > Rº dada por 


n—l 
gm) = +55 xy, 
i=1 
Como À é aberto existe à, >0 tal que B,1 (p) c A. Mas g é contínua e g(0) = p, logo 
, e obtemos f: B;(0) > A. Ainda, 
fé de classe Cº e, como DF (0) e X(f (0)) são linearmente independentes, f é uma seção 


existe ô >0 tal que g(B;(0)) c B,, (p). Fazemos f= dla; o 
transversal local de X satisfazendo f (0) = p. Temos ainda que > é uma seção transversal 
de X em p. 


Teorema 2.5.3 (Fluxo Tubular) Seja p um ponto regular de X: A > Rº de classe C* 
ef: A> > uma seção transversal local de X de classe C* com f (0) = p. Então existe uma 
vizinhança V de p em A e um difeomortismo h: V > (-e, e) x B de classe C*%, ondee >0 e 
B é uma bola aberta em R"! de centro na origem 0 = f-'(p) tal que: 

LL hHEnV)=(0)xB; 
Il h é uma C'-conjugação entre XI, e o campo constante Y: (-e, e) x B> Rº, 


Y=(1,0,0,:-:,0)€ R”. 


Teorema 2.5.4 Sejam > uma seção transversal de um campo X de classe C', definida 
no aberto Ade R", e pe >. Então existem uma vizinhança V de p, e= e(p) >0 e uma função 
T:V>(-e ejtalquerT(Vn>)=0€e: 


a) Para qualquer qe V, a curva integral q(:, q) de XI, é definida e biunívoca 
emJ,=(T(g)-e T(9) + 6); 


b) n(g) := p(T (a), q) é o único ponto onde q (:, ql, intercepta a 2; 


c)n: V> é de classe C*e Dn(g) :R" > T, 92 é sobrejetiva para qualquer q € 


V. Mais ainda, Dn, -V=0 se, e somente se, v= ax(g) para algum ae R 
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2.6 ESTABILIDADE DE LIAPUNOV 


Nesta seção, trataremos do Critério de Liapunov para estabilidade de um ponto 
singular de um sistema autônomo. Mais geralmente, consideremos o campo dado pelo 


sistema 


onde f: A => Rº é de classe C!. 

Definição 2.6.1 Um ponto singular p € A do campo dado por (2.2) é dito estável 
quando para toda vizinhança U de p, existe uma vizinhança U, de p tal que toda solução q 
de (2.2) com qp(0) E U,, tem-se que p(i) E U para todo t> 0. Se além disso dim o (d =p, 
então p é dito assintoticamente estável. 

Seja V: A > R uma função diferenciável. Se y é solução de (2.2) definida no 
intervalo |, então 
(V(g())=Viv(d) pd dt 

=VV(o(d) fold), tel, 


dt 


onde - denota o produto escalar de R". Denotando por x = q (0) temos 


d 
ZE Vigo (o =VVO) fl). 


Ponhamos, para cada x€E À, 
V= VV (9 - f(. 


Definição 2.6.2 Seja p € A um ponto singular de (2.2). Uma função de Liapunov para 
p é uma função V: U> R, de classe C', definida em um aberto U contendo p, satisfazendo 
as seguintes condições: 


L V(p=0eV(M)>0;Vx= p; 
H V<0emu. 


A função de Liapunov V diz-se estrita quando 
ll. V<0emU-(p). 


O critério de Liapunov para o sistema (2.2) é dado pelo teorema seguinte: 


Teorema 2.6.3 Seja p um ponto singular de (2.2). Se existe uma função de Liapunov 
para p, então p é estável. Se a função for estrita, então p é assintoticamente estável. 
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2.7 TEOREMA DE POINCARÉ-BENDIXSON 


Definição 2.7.1 Sejam A c R" um aberto, e X: A > R” um campo de vetores de 
classe C&, k=1,2,...,º ouk=w. Seja 1, =(w (0), w (p)) O intervalo máximo de pe A. Se 
w (p) =ºº, define-se o conjunto 

w(p) =tge A;3(t) comi > ve p(t, p) > q, quando n > 00). 

Analogamente, se w (p) = -º, define-se o conjunto 

a(p) =tge A, a(t) comt, > -v e p(t, p) > q, quando n > 00). 

Os conjuntos w(p) e a(p) são chamados, respectivamente, de conjunto w-limite e 
conjunto a-limite de p. 

Exemplo 2.7.2 Seja X: R? > Rº dado por X(x, y) = (x, —y). Assim, q (t p) = (ce, 
ce). Consequentemente, 


Se p = (0, 0), então a(p) = w(p) = (0, 0)». 
Sepet(x,0); x 0), então a(p) =((0, 0); e w(p) = 6. 
Se pe (O, y); y * 03, então a(p) = d e w(p) =((0, 0)). 
Sepetf(x, y); xy * 0), então a(p) = w(p) = 6. 

Veja Figura 1.2. 


Observação 2.7.3 Sejam q, a curva integral do campo X pelo ponto pe y, a 
curva inte- gral do campo —-X pelo ponto p, então P -) = vd, vt e | Segue daí que o 
w-limite de p com relação ao campo X é igual ao a-limite de p com relação ao campo —X, 
e reciprocamente. Por este motivo, para estudarmos as propriedades gerais dos conjuntos 
a-limite e w-limite de órbitas é suficiente nos restringirmos ao estudos do conjunto w-limite. 
Apresentamos a seguir propriedades gerais do conjunto w-limite. 


Teorema 2.7.4 Sejam X: A > R" um campo de vetores de classe C', evy(p) = (op (t, 
p); t> 0) (respectivamente y(p) = (y(t, p): t<0)) a semi-órbita positiva (respectivamente 
a semi-órbita negativa) do campo X pelo ponto p. Se y*(p) (respectivamente y(p)) está 
contida num subconjunto compacto K c A, então: 
| w(p) + 6 (respectivamente a(p)); 


HI. o(p) é compacto (respectivamente a(p)); 
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Figura 2.2 — Retrato de fase do sistema do Exemplo 2.7.2. 


Fonte: Elaborado pelo autor. 


HW. o(p) é invariante por X (respectivamente a(p)), isto é, se q E w(p), então a 
curva integral de X por q está contido em w(p); 


IV. w(p) é conexo (respectivamente a(p)); 
Vo dt p), w(p)) > O quando t > +º. 

Observação 2.7.5 Se y, é a órbita de X pelo ponto pe qe y,, então w(p) = w(q). 

Com efeito, se q E A existe ce Rtal que q (t p) = q (t+c, q). Analogamente, a(p) = 
a(g). A observação anterior permite que definamos o conjunto w-limite w(I”) de uma órbita 
F como o conjunto w-limite de qualquer um dos pontos de T. 

Nosso próximo objetivo é apresentar o Teorema de Poincaré-Bendixson. Trata-se 
de um dos principais resultados da Teoria Qualitativa de EDO's no plano, pois classifica os 
conjuntos w-limites de campos planares. Cabe ressaltar que o resultado não é válido para 
campos de dimensão > 3. Lembremos que a semi-órbita positiva por p, denotada por y*, 
segundo o Teorema 2.7.4, é definida por y; = (q (t, p); t > 0). 

Definição 2.7.6 Dado um campo de vetores X: A > Rº, onde A é um aberto de Rº, 
entendemos por ciclo limite uma órbita periódica de X que seja isolada no conjunto de todas 
as órbitas periódicas de X, ou seja, uma órbita periódica y é um ciclo limite se existe uma 
vizinhança V de tal forma que y seja a única órbita periódica contida em V. 
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Teorema 2.7.7 (Poincaré-Bendixson) Sejam X: A > R2, campo de classe C', 
definido no aberio A c R2e pe A tal que V esteja contida num compacto K c A. Suponha 
que o campo X possua um número finito de singularidades em w(p). Têm-se as seguintes 
alternativas: 

| | sew(p) contém somente pontos regulares, então w(p) é uma órbita periódica; 


H. se w(p) contém pontos regulares e singulares, então w(p) consiste de 
um conjunto de órbitas, cada uma das quais tende a um desses pontos 
singulares quando t > +90; 


ll. se w(p) não contém pontos regulares, então w(p) é um ponto singular. 


Corolário 2.7.8 Seja X um campo de vetores de classe C! num conjunto aberto A 
c R2. Sey é uma órbita fechada de X tal que Inty c A então existe um ponto singular de X 
contido em Inty. 

A demonstração do Teorema de Poincaré-Bendixson é feita em uma sequência de 
lemas. Destacamos um deles, que será útil no decorrer do texto. 

Lema 2.7.9 Seja X um campo de vetores de classe C! e um conjunto aberto A c Rê, 
e > uma seção transversal a X contida em À. Se y é uma órbitade Xe pe X N y, então, vi= 
to(t p); t>0) intercepta > numa sequência monótona p,, D,, --., Py -- 


2.8 EXISTÊNCIA DE ÓRBITA PERIÓDICA PARA A EQUAÇÃO DE LIÉNARD 

Utilizaremos o Teorema de Poincaré-Bendixson para demonstrar a existência de 
uma órbita periódica para a equação de Liénard, uma das mais estudadas equações 
diferenciais planares. 

No Capítulo 3, estudaremos a equação de Van de Pol, um caso particular da equa- 
ção de Liénard. Nessa oportunidade, desenvolveremos uma ferramenta que permite che- 
gar de modo mais conciso a resultados mais específicos sobre a equação de Van de Pol. 

Nesta seção seguiremos os resultados de J. SOTOMAYOR[09] e L. PERKO[10]. 


Teorema 2.8.1 [Equação de Liénard] Seja g: R > R uma função de classe C! tal 
que G(u) =Jug(s) ds é ímparem u, G(u) > seu > o, existeB>0O tal queseu>B, Gé 
crescente, e existea >0 tal que G(u) <0 se0O<u<a. Então, a equação de segunda ordem 


u"+ g(luu'+u=0 (2.3) 


que é equivalente ao sistema 
E =v—G(u), 


v=—u. (2.4) 


admite uma solução periódica. 
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Demonstração: Mosiraremos que, dado v, suficientemente grande, o conjunto 


w(A), onde A = (0, v,), está contido em uma região limitada, e através do Teorema de 


Poincaré- Bendixson, é uma órbita fechada. Mas antes, observemos algumas propriedades 


do sistema (2.4). 


1. 
2 


A única singularidade de (2.4) é 0 = (0, 0). 


Para a solução t » (u(i), v(?)) de (2.4) temos que: u é crescente onde v(f) > 
G(u(t)) e decrescente onde v(?) < G(u(t)). Também v é decrescente se u(f) >0 e 
crescente se u(t) < 0. Além disso, o campo (v — G(u), —u) é horizontal no eixo v 
e vertical na curva v= G(u). 


Do fato de G ser ímpar, decorre que as soluções t » (u(?), v(?)) de (2.4) são in- 
variantes por reflexões, isto é, t » (u(?), v(?)) é solução de (2.4) se, e somente 
se, tr (-=u(i), —v(?)) também o for. Em particular, conhecendo o arco ABCD 
de uma órbita de (2.4), como mostra a Figura 1.3, sua reflexão com respeito 
a origem também é arco de uma órbita. Se A= (0, v), D=(0,-v)e v,< 
v, então a semi- órbita positiva que passa por A é limitada e portanto está 
contida na região limitada pela curva de Jordan J formada pelo arco ABECD, 
por sua reflexão com respeito à origem e pelos segmentos do eixo v que ligam 
os extremos desses arcos. 


Aqui, A=(0, v), B=(u, v), E=(u, v)), C=(u, V;), D=(0, v,) e y = ABECD é um 
arco da trajetória de (2.4). 


De fato, vamos mostrar que qualquer solução de (2.4) saindo do ponto A = (0, v,), 


com v, >? poisv <-us-B > v() >-Bt+ [B?, tem uma órbita tal como mostra a Figura 


1.3. Vamos verificar que as soluções se aproximam da origem, para isso observemos as 


curvas onde tt =0e )=0. 


Definamos: 


V+=t(0,v): v>0% 
V-=((0, o v<0), 
g=u, v):v= G(u), u>0), 
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Figura 2.3 — Retrato de fase do sistema (2.4). 


Fonte: Elaborado pelo autor. 


g=(u, v):v=G(u), u<o0Oy, 

e consideremos Z a região entre V+e gt, E a região entre gre V-, K a região entre 
V-eg,e0aregiãoentreg e V+. 

Sejam (u, V)) E V+e (u(i), v(1)) solução de (2.4) passando por (u,, v.). Logo, u'(0) >0 
de modo que a solução entra em Ze u, = u(t,) >0, para t, >0 pequeno, enquanto a solução 
permanece na região g. Assim, se 1 >0 então u( > u.. 

Seja a região ((u, y:uzu, v<v, (u, VE G com, U <-u<o0, então a solução só 
pode sair de q atravessando a curva g* em algum momento t, >0. 

Seja (u,, V,) E g' a solução no tempo t,. A trajetória entra em € para algum tempo t, 
>. 

Além disso, ao longo de todo g*, v <0; então uma vez que a trajetória sai de (u,, 
V;), ela não pode retornar a g*. Assim, ao longo da trajetória, há um limite superior em ii, ti 
s-a<0parat>t, e u(?) = u, - a(t- t,) enquanto permanece na região €. 

Então, a trajetória pode deixar & apenas passando por V- ou tornando-se ilimitada. 
Por tw <-a<o0, ut) se torna zero quando t = t, + u,/a. No entanto, nesse intervalo de 
tempo, v=-u>-u, eentão Vl) >v,-u, (ua). 

Mas dado que a priori, v(t) é limitada nesse intervalo, a solução deve atravessar 
V- saindo de é. 

Pela simetria da equação, podemos fazer a mesma argumentação, segundo ajustes 
apropriados, para a região simétrica ao eixo v. 
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Agora vamos mosirar que se v, é suficientemente grande temos que v, < v, e 
portanto o conjunto O(A) estará contido em uma região limitada, e então, pelo Teorema de 
Poincaré-Bendixon, O(A) será uma órbita fechada já que O(A) = (0, 0) pois, (0, 0) é uma 
fonte de (2.4), como mostraremos abaixo. 

Seja a função 

R(u, v) = - (u + v). 


4 


Para uma solução u = u(?), v = v(?) de (2.4) temos 
= (ul, VC) = u(v — G(u)) + v(-u) = —u() G(u(t)). (2.5) 


Assim temos, com A=(0, v)e D=(0, v,) 


1 232 2 = Es 
03 —v8) = RD) — RA) , oo R 


E Ia E A Res Lo ar. 


Pelo Lema 2.8.2, segue 


5 (vá 9 =| [ / ] Rn dia A Glu)do 


Pelo Lema 2.8.3, segue que as duas primeiras integrais tendem monotonicamente 


a zero quando v, > 00, 
Ao ca do segmento BEC os limites de integrção correspondentes a v estão fixos 


emv=v ev=v, eparacada ve [v,, v,) temos então 


1a 


va 
/ Glujdo = | G(u)dv 
BEC v] 


ú(vo) = / G(u)do 


E! 


Notemos que 


satisfaz a 


ú(vo) “eu 


vi 


Logo, para € > 0 suficientemente pequeno , 


vi vi+e 
-o(m) = [ Glujdo > [ G(u)dv 


a 3—€ 
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> G(e) ih “do = G(e)[vy — vg — 2€] > G(o)[vy — 2]. 


va—e 
Como v, > quando v, > 9º, isto prova que W(v,) > -o se v, > o, Portanto vã < vê, 
se v, é grande. 
Assim, por (2.5), se 0 <lul <a, R (u(f), V()) >0, logo -R é definida negativa. 
Então, pelo Teorema de Estabilidade de Liapunov, (0, 0) é solução assintoticamente 
estável de (2.4) para t > -o, e portanto (0, 0) é uma fonte de (2.4), isto é, (0,0) é o a-limite 
de todo ponto numa vizinhança de (0, 0), terminando aqui a prova de que (0, 0) é uma fonte. 
Lema 2.8.2 Seja R como na demonstração do Teorema 2.8.1. Então, 


[Lo* Lo di a im / + ] a du + no G(u)do 


Demonstração: Em AB consideremos a parametrização u > (u, v(u)) e teremos 


Análogo para CD. 


Já para BEC consideremos a parametrização v > (u(v), V)), V, < V=< v, e temos 


1 du dt 
(u(v),v) = um +u=u ET do +u 


= u(v— el) E) +v 


= G(u(v)). 


Assim, 


vz 
-[ di = ar = [ G(u(v))do. 
BEC CEB vi 


Lema 2.8.3 Sob as mesmas hipóteses do Lema 2.8.2 para u, v e G(u), temos que 


INE 


tendem monotonicamente a zero quando v, — 00. 

Observemos que a igualdade que decorre do Lema 2.8.2 para as duas integrais tem 
o mesmo intervalo de integração e a única ressalva é que uma integral é para a curva de 
cima, enquanto a outra é a integral pela mesma curva por baixo. 
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Demonstração: Sejam M=maxtu+ G(u):0<u<fB)e N=maxluG(u)l:0<u<B)e 
tomando v, suficientemente grande de modo que v(8) >M, então no intervalo considerado 


temos: 
dv —u 
dE. se vlw) > vn —“ 
du  v—G(u) PEA MA 
Logo, 
2 uG(u) 2 JuG(u)| Po JuG(u)] 
— ++ du! < — ++ du < — > ww >——w+w+ du. 
/ vu) — Gu) ] a v(u) — G(u) < [ vo — (u+G(u)) 
Assim, 


o 
O<I/ — (o mc agtlsN [ 4 
v(u) — vo — 


quando v, > 0, 
Portanto, temos que [,, dR > O quando v, > o. Prova análoga para [.,dR > 0 
quando v, > ºº, pois v, > o implica v, — 00, 


2.9 INTEGRAIS PRIMEIRAS E SISTEMAS HAMILTONIANOS 


Definição 2.9.1 Sejam X um campo de vetores de classe C' definido em um aberto 
AdeRºeH:A> R uma aplicação de classe C!. Dizemos que H é uma integral primeira de 
X em A se H permanece constante ao longo das órbitas de X contidas em A e, além disso, 
H não é constante em todo subconjunto aberto de A. 

Desta forma, Hétal que (H º qp)(t)=0,te |, para toda solução q = q (t) do campo 
X definida no intervalo | c R. Em particular, cada solução do campo permanece em uma 
única curva de nível H. = ((x, y) E A; H(x, y) = c). Assim, as curvas de nível H. em A são 
suficientes para descrever o retrato de fase de X em A. 

Definição 2.9.2 Sejam A um subconjunto aberto doR? e H: A > R uma aplicação de 
classe C2. Um sistema da forma 


0H 9H 
=>, y=->>, (9) € A (2.6) 


é chamado um sistema hamiltoniano. 

A aplicação H é cnamada uma função Hamiltoniana do campo dado em (2.6). 

Note que H é constante ao longo de cada solução de (2.6). De fato, dada uma 
solução q (?) = (x(D, Y()) definida em um intervalo | c R de (2.6) temos: 
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Consequentemente, para descrever o retrato de fase de um sistema Hamiltoniano 
em Rº basta descrever as curvas de nível H(x, = c, ce R. 

O exemplo que segue trata de um sistema que descreve o movimento de um pêndulo 
não-amortecido, por ser equivalente à equação de segunda ordem % = -kx. Observe que a 
função hamiltoniana nesse caso é a energia total do sistema. 

Exemplo 2.9.3 Considere o sistema 

t=y y=-—kax 


e a função H dada por 
1 
H(x,y) = a + ka?), 


onde k >0 é uma constante. Segue que 


0H 9H 
da (ou) = he e go(m 8) = (0,9) E Rº. 
Assim, o sistema acima pode ser reescrito da seguinte forma 
a = O 0H 
=", U== 
Oy dx 


ou seja, é Hamiltoniano. 

Em seguida, mostraremos que os fluxos de sistemas Hamiltonianos preservam 
área, no seguinte sentido: 

Definição 2.9.4 Um difeomorfismo h: A > A, onde A é um subconjunto aberto do 
R”, preserva área quando, para qualquer aberto A, c A de área finita, os conjuntos h(A,) e 
A, têm a mesma área. 

Teorema 2.9.5 Um difeomorfismo h: A> A, de classe C!, com Ac Rº preserva área 
se, e somente se, o determinante de Dh tem valor absoluto constante e igual a 1. 

Definição 2.9.6 Seja X: A > R2, um campo de vetores de classe C' definido num 
aberto Ac Rº cujo fluxo q está definido para todo (t, x) e R x A. Dizemos que o fluxo do 
campo X preserva área se, para cada te R e todo aberto À, & A de área finita temos 


sendo q, a aplicação dada por x » q (t, x). 
Proposição 2.9.7 (Fórmula de Liouville-Ostragradski). Se X: A => R? é um campo de 
classe C!, no aberto Ac Rºe q oseu fluxo, então: 


detD(pua) = exp ( / | dio(X(pslo)Dás) 


com xe a. 
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Teorema 2.9.8 (Teorema de Liouville). Se um campo de vetores planar de classe C! 
tem divergente nulo, então seu fluxo preserva área. 

Notemos que o Teorema 2.9.8 nos dá algumas relações entre as trajetórias de um 
campo de vetores planar, pois considerando um feixe de trajetórias com condições iniciais 
em um aberto A c R2, como esse fluxo preserva área não é possível que todas essas 
trajetórias se afastem (respectivamente, aproximem) umas das outras, já que isso causaria 
uma expansão (respectivamente, compressão) da área inicial. 

Proposição 2.9.9 Fluxos associados a sistemas Hamiltonianos preservam área. 
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CRITÉRIOS DE EXISTÊNCIA E NÃO EXISTÊNCIA DE CICLOS 
LIMITES EM SISTEMAS PLANARES 


O conceito de ciclo limite para campos vetoriais planares, introduzido por Poincaré, 
modela oscilações não lineares em processos físicos. No final da década de 1920, Van 
der Pol, Liénard e Andronov, estudando oscilações não lineares de fenômenos elétricos, 
analisaram certas equações diferenciais ordinárias de segunda ordem e verifica- ram 
a ocorrência de ciclos limites. Após tal verificação, matemáticos e físicos estudaram 
extensivamente a existência e unicidade, entre outras propriedades, destes ciclos limites. 
Este capítulo está baseado no Capítulo 7 de F. DUMORTIER, J. LLIBRE and J.C. ARTÉS 
[01]. 


3.1 RESULTADOS BÁSICOS 


Consideramos o sistema de equações diferenciais 


o dx od 
i=D =P(x,y), = — U,y), a 


onde xe y são variáveis reais, Pe Q são funções C'! definidas em um aberto A de 
R2. Assim, estamos sob as condições do Teorema de Existência e Unicidade de Picard. 

Classificamos a seguir os diversos tipos de órbitas periódicas de campos planares. 
Seja y uma órbita periódica do sistema (3.1). Se para alguma vizinhança externa (interna) 
da órbita periódica y não existem outras órbitas periódicas de (3.1), então y é chamado de 
ciclo limite externo (interno). Uma órbita periódica que seja ciclo limite externo e interno é 
chamada de ciclo limite. 

Assim, um ciclo limite é uma órbita periódica isolada de (3.1). 

Se existe uma vizinhança externa (interna) da órbita periódica y formada pela união 
de órbitas periódicas, então y é cnamada de órbita periódica de tipo centro-externa (centro- 
interna). 

Se em qualquer vizinhança externa (interna) de y existem órbitas periódicas e não 
periódicas, diferentes de y, então y é chamado de um ciclo externo (interno) indefinido. Se 
uma órbita periódica y for ciclo externo indefinido ou ciclo interno indefinido, chamamos y 
de um ciclo indefinido. 

Proposição 3.1.1 Se y é uma órbita periódica do sistema planar (3.1) então existe 
uma vizinhança U suficientemente pequena dey satisfazendo as três seguintes afirmações: 


* (a) Unão contém pontos singulares; 
*— (b) para qualquer ponto p E y existe um segmento transversal > á de X em p; 


* (oc) toda órbita periódica contida em U intercepta >, transversalmente em um 
único ponto. Qualquer órbita não periódica em U intercepta >, em um número 
infinito de pontos, que estão no mesmo lado de y. 
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Demonstração: 


a. Seja p € y qualquer, então X(p) + 0, pois y é órbita periódica. Pela continui- 
dade, existe aberto U, tal que q E U, X(q) + O. Basta tomarmos U = at 
assim U é aberto que não contém pontos singulares. 


b. O campo X é de classe C', U, é aberto, logo pela Definição 1.5.1 temos que 
para cada p E y, existe um segmento transversal 2, 


c. Sejam 2, uma seção transversal de Xem p com pe y,e Va vizinhança aberta 
de E, dada pelo Teorema do Fluxo Tubular. Então, toda órbita de X que passa 
por Vintercepta > o transversalmente. 


Para cada te R definamos 
W=to(t x;xeV), 

que é um conjunto aberto pela Observação 2.3.8. 

Seja T>0 período de y, e definamos 

w=Jw 
te[0,7] 

que é um aberto do Rº2. 

Diminuindo U se necessário, podemos supor que Wc U. Dado gq E U, temos da 
definição de W que q = q (t, x), para algum te [0, 7], e algum xe 2,. Logo, a órbita que 
passa por gintercepta ê, em um ponto, a saber, o ponto q (-t, q) = x. Como toda órbita que 
passa por Vintercepta Ê, transversalmente, o aberto de g intercepta a transversalmente. 

Se Y, não é uma órbita periódica, pelo Lema 2.7.9, v; intercepta 2, segundo uma 
sequência monótona de pontos. Se dois desses pontos coincidem, então Y, Seria perió- 
dica. Se dois desses pontos estivessem em lados distintos de y então, por continuidade, 
algum ponto de Y interceptaria y, contradizendo o Teorema de Picard. 


3.2 A TRANSFORMAÇÃO DE POINCARÉ 


A ferramenta básica para o estudo de órbitas periódicas é a Transformação de 
Poincaré associada a uma órbita periódica y, a qual é um difeomorfismo 7 que definiremos 
a seguir. Esta transformação descreve o comportamento das órbitas do campo numa 
vizinhança de y. 

Apresentaremos a definição da Transformação de Poincaré para campos em Rº, e 
depois a utilizaremos para estudar ciclos limites de sistemas planares. 

Seja então y = (y(t p); O <t< Tuma órbita periódica de período T, > 0 de um 
campo X de classe C*, k > 1, definido em Vc R". Seja > uma seção transversal de X em p. 
Em virtude da continuidade do fluxo q de X, para todo ponto qg E > próximo de p, a trajetória 
t+ q(t, q) permanece próxima a y, com tno intervalo compacto [0, 27]. Define-se 7 (q) 
como o primeiro ponto onde y volta a interceptar >. Seja >, o domínio de 7. Naturalmente 
peZ,e 7(p)=p, Veja Figura 2.1. 
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A seguir, demonsiraremos que 7 : >, > > é um difeomorfismo de classe C* sobre 
sua imagem >,. Vamos usar o Teorema do Fluxo Tubular. Como q (T,, p) = p, existe uma 
vizinhança 2, de pem X tal que y(T,, q) E Vparatodo qe 2,. Sejan: V-> > a aplicação 
definida já definida, e pomos 7 :>,>2, m(g)=n(p(T, q). 

Outra expressão para 77 é m(q)= q(T,+ T(P(T, Q)). 9), onde T: V=> R é o tempo 
T (x) que leva a órbita por xem V para interceptar >. Do Teorema da Função Implícita, T é 
de classe C*. 

Desta expressão resulta que 7 é da mesma classe de diferenciabilidade que X. A 
inversa 77!:2, >2, de 7 é definida tomando-se a campo -X. Fica provado que 7 é um 
difeomorfismo C*, 

Muitas propriedades do retrato de fase de X perto de y se refletem em 71 e recipro- 
camente. Por exemplo, as órbitas periódicas de X vizinhas de y correspondem aos pontos 
fixos de 7, que são pontos g€ 2, para as quais 71 (q) = q. O comportamento assintótico das 
órbitas de X perto de y também é descrito por 7. Assim, é possível demonstrar que lim, 
m."(q) = p implica lim,.. (gv (t q), y) = 0, onde dp (t q), y) = inlgp (t q)-n, rey). 

Definição 3.2.1 Com a notação acima, a órbita fechada y é um atrator periódico 
(ou então y diz-se orbitalmente estável) quando lim, .. d(p (t q), y) = O para todo q numa 
vizinhança de y. 

O próximo resultado mostra que existem apenas três tipos de ciclos limites em sis- 
temas planares. 

Proposição 3.2.2 Sejam X: V> Rº um campo vetorial de classe C*, y um ciclo limite 
de Xe V, a vizinhança de y tal que y é a única órbita periódica de X que intercepta V,. Então 
exitem apenas os seguintes tipos de ciclos limites: 


Figura 3.1 — Seção transversal > e transformação de Poincaré. 


> 


Fonte: Elaborado pelo autor. 
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(a) estável, quando lim, (vp (t q), )=0 paratodoge V; 
(b) instável, quando lim, .. d(p(t q), y) =0 paratodoge V; 


(c) semi-estável, quando lim, ., d(y (t, q), ) = 0 para todo ge Vn Exty, elim,. .. 
d(p(t q), ) =0 para todo ge Vn Inty, ouo contrário. 

Veja Figura 2.2. 

Demonstração: Seja p € y qualquer. Diminuindo a vizinhança V, se necessário, 
podemos supor pela Proposição 2.1.1, que ela não contém singularidades e que > é uma 
seção transversal a Xem p. Seja 7 : >, > > a transformação de Poincaré. Suponhamos 
que > esteja ordenado, sendo o sentido positivo de Exty para Inty. Dado qe 2, N Exty, 
temos que 77 (9) + q pois, caso contrário, a órbita que contém q seria periódica. Então 71 (q) 
>qou (9) <q. 

Suponhamos 7 (q) > q. Seja agora a região A limitada por y, pelo arco de trajetória 


= 


ar(q)e pelo segmento gr(q) E X,. A região A é compacta, logo pelo item v) do Teorema 


2.7.4 temos que d(q (t, p), w(p)) > O quando t > +00, 

Se 7 (q) < q, considerando o campo —X, fica provado que dp (t p), w(p)) > 0 
quando t > +90. 

As mesmas considerações podem ser feitas em Inty. Combinando todas as possi- 
bilidades podemos provar a proposição. 


Figura 3.2 — Classificação de ciclos limites. 


O S 


Estável , 
Instável 


O S 


Semi-estável 


Fonte: Elaborado pelo autor. 
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Na notação da Proposição 3.2.2-(c), quando ocorre lim, ., d(y(t, q), y) = O para 
todo qe Vn Exty, dizemos que y é um ciclo limite externamente estável. Análogo para 
internamente estável, externamente instável e internamente instável. 

Observação 3.2.3 Com as notações da Proposição 3.2.2, segue que y é um ciclo 
limite se, e somente se, p é um ponto fixo isolado de rt. Além disso, 


*  yéestávelse, e somente se, Ir (x) — pl <Ix— pl para todo x + p suficientemente 
próximo de p; 


*- yéinstável se, e somente se, lt (x) — pl >Ix-— pl para todo x  p suficientemete 
próximo de p; 


*- y é semi-estável se, e somente se, Ir (x) — pl <Ix— pl para todo x E Xn Exty 
suficientemente próximo de p el7 (x)- pl >Ix- pl para todo x € 2n Inty suficien- 
temente próximo de p, ou o contrário. 

Proposição 3.2.4 Qualquer órbita periódica y é, ou um ciclo limite externo (interno), 
ou uma órbita periódica do tipo centro-externa (interna) ou ainda um ciclo externo indefinido 
(interno). 

Demonstração: Consideraremos a transformação de Poincaré 7 : >, > >. Então 
existem 3 casos. 

Primeiro caso: existe e >0 tal que 71 (x) =x, para todo x E (p, p + €). Esse caso segue 
imediatamente da Proposição 3.2.2. 

Segundo caso: existe e >0 tal que 7 (x) = x para qualquer x E (p, p + €), é uma órbita 
periódica do tipo centro-externa (interna), pois existe uma vizinhança externa (interna) da 
órbita periódica formada pela união de órbitas periódicas. 

Terceiro caso: para todo €>0, existem x, ye (p,p+ejonde T(x) =xe T(J) = y. 

Nesse caso, da definição da transformação de Poincaré, a órbita de x é periódica e 
a órbita de y não. 


Figura 3.3 — Órbita do tipo centro-externa e ciclo externo indefinido. 


Fonte: Elaborado pelo autor. 
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A Figura 2.3 mostra os casos 2 e 3 da Proposição 3.2.4. 

Observação 3.2.5 Como 71 é diferenciável, se 77'(p) < 1 podemos concluir que y é 
estável e se 7'(p) >1 que y é instável. 

Demonstramos a seguir que, para os sistemas analíticos, ocorrem uma restrição na 
classificação das órbitas periódicas. 

Proposição 3.2.6 Se Pe Q são funções analíticas, então o sistema (3.1) não tem 
ciclos indefinidos. 

Demonstração: Por Pe Q serem funções analíticas, temos que a aplicação de 
Poincaré também é analítica. Então pelo princípio dos zeros isolados das funções analíticas, 
temos que o sistema (3.1) não ocorre o terceiro caso da demonstração da Proposição 3.2.4. 

Definição 3.2.7 O divergente do campo X = (P Q) é o traço da matriz jacobiana de 
X, ou seja, 


divd(X(x, y)) = trDX(x, y) = ua (x, y) + dq 


dr D) Oy (x, y). 


O teorema a seguir nos dá condições sobre a estabilidade de ciclos limites de sis- 
temas planares analisando o divergente do campo. 

Teorema 3.2.8 Sejam X=(P Q): V> Rº um campo vetorial de classe C!, y uma 
órbita periódica de X de período Te 7:22, > > a transformação de Poincaré associada a 
uma seção transversal > de Xem pe y. Então: 


"o = es ( [ divexto(opa) e (( (Soto) + Soo) a). 


Em particular, seJ7 divX(y())dt < O então y é estável e seJ7 divX(y(i))dt > 0, então 
Y é instável. 

Demonstração: Para cada t, suponhamos A(t) = DX(v(1)). Seja ó (1) a matriz funda- 
mental de x = A()x, com &(0) = 1: Pela fórmula de Liouville, temos 


det (7) = exp [/ dioX (tod 


Vamos mostrar que 7'(p) = detó(T). Seja q o fluxo de X. Temos que d(T) = D,qp(T, 
p), e notemos primeiro que D, q (T, p).X(p) = X(p). De fato, como D q(t, pJl,, = X(p), vem 


AD lt) lco = LolT + tipo 


p(t, p)lt=o — X(p). 


Doy(T,p).X(p) = Prá 


d 
=” 


Por outro lado, se g : (-€, e) > X é uma parametrização de > tal que g(0) = p, O 
conjunto B = (X(p), g' (0); é uma base de R?, Por definição, r(g(s)) = q(T+ T(W(T g(s))), 
g(s)), donde 
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a 


qe” º 9(S)|s=o = Di(T,pja + Do(T, p).g'(0) 
= aX(p) + Dop(T,p).g(0), 


onde a é a derivada de T (q (T, g(s))) em s = 0. Portanto, a matriz de D,q(T, p) na 


m'(p).g'(0) 


base B é 


e assim obtemos que detó(T) = 7'(p). 


3.3 EXISTÊNCIA E NÃO-EXISTÊNCIA DE CICLOS LIMITES EM CAMPOS 
PLANARES 

Nesta seção apresentamos resultados clássicos relacionados à existência e não- 
existência de ciclos limites em campos vetoriais planares. 

Quando nos referirmos a uma região anular, estamos considerando como sendo 
a região delimitada por duas curvas simples fechadas do plano, com interseção igual ao 
vazio, e que uma esteja contida no interior da outra. 

Proposição 3.3.1 Se y, e y, são órbitas periódicas que formam a fronteira de uma 
região anular A, e não existem pontos singulares e nem outras órbitas periódicas em A, 
então a órbita y, em Aé o conjunto w-limite (a-limite) e y, é o conjunto a-limite (w-limite) de 
todos os pontos de A. 

Demonstração: Sejam y, e y, órbitas periódicas que formam a fronteira de uma 
região anular A; obviamente A é compacto. Seja p E 4. Então, pelo Teorema de Picard, Y 
c A. Logo w(p) db. 

Pelo Teorema de Poincaré-Bendixson, w(p) é órbita periódica e suponhamos w(p) = 
Y, Então, pelo Lema 2.7.9, teremos que a(p) = y,. Resta mostrar que Vg € 4, vale que w(q) 
=y,e a(g) = y,. Se w(9) = y, então temos uma contradição com a continuidade do fluxo. 


Teorema 3.3.2 (Região Anular de Poincaré) Seja A uma região anular que não 
contém pontos singulares, dA não contém órbitas periódicas e que cada órbita que cruza a 
fronteira de A move-se do exterior para o interior. Então existe pelo menos um ciclo limite 
externamente estável em A e existe também um ciclo limite internamente estável em A, 
sendo possível que ambos ciclos limites coincidam, de modo a ser um único ciclo limite 
estável. 

Demonstração: A ideia para a prova consiste em utilizar o Lema de Zorn. 
Consideremos a seguinte família de conjuntos 

F=(y:y órbita periódica contida em A). 

Pelo Teorema de Poincaré-Bendixson temos que T' + 6, e estabelecendo a seguinte 

relação 
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Y4,SYO Tnty E mta 
facilmente vemos que < é uma relação de ordem (reflexiva, antissimétrica e transitiva) 
de T. Seja agora S um conjunto totalmente ordenado de T', devemos mostrar que S tem 
uma cota superior. Seja 
d= N [Inty;N A]. 


NES 


Notemos que ô + & pois cada Inty;N A é compacto e a família (Tnty, NA, y,€ S) 
tem a propriedade da interseção finita. 

Seja q € 0. Temos então que w(q) é uma órbita periódica pelo Teorema de Poincaré- 
Bendixson e w(9) E O (pois O é invariante e não contém pontos singulares), e assim, O(g) 
é cota superior de S. 

Pelo Lema de Zorn, [ contém um elemento maximal M, comy<M> IntM c 
Inty, Vy ET. 

Então, M é órbita periódica internamente estável. 

A outra parte da demonstração é inteiramente análoga, fazendo os devidos ajustes 
dos conjuntos. 


Proposição 3.3.3 (Princípio da Simetria) Suponhamos que no sistema (3.1) ocorre 
P(x, y) = P(-x, y), Mx, y) = -Q(x, y), V(x, y) E Rº, 


e a origem é o único ponto singular no eixo y. Se a órbita y começa sobre o eixo 
positivo y e retorna ao eixo y negativo, então y é uma órbita periódica. Se todas as trajetórias 
próximas à origem possuem essas propriedades, então a origem é um centro. 

Observe que a hipótese de que a órbita comece e retorne ao eixo y é essencial, 
como mostra a Figura 2.4. 

Demonstração: Digamos que y seja a imagem da trajetória f(?) = (x(d, y(d) para t 
E |, sendo [intervalo máximo. 

Suponha que f(0) = (0, y)) e f(t) = (0, —y,) para algum £, >0, sendo y,, y, >0. 

É suficiente mostrar que F(-t) = f(t). Para isso, provamos inicialmente que —t € 1. 
Definamos 


o =-x(-), com te -I=(sE R;-sE | 
VD = y(-D. (3.2) 


Então, por hipótese, verifica-se que: 
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ut) = P(u(t), v(t)) 

(8) = Q(u(t), (8), 

u(0) = 0, (3.3) 
v(O) = yo 


Pelo Teorema de Picard, temos u(t) = x(?) e v(?) = y(t), para qualquer te —|, logo x(-t) 
=-x(?), y(?) = y(-?), para qualquer te —I. 

Assim, multiplicando a primeira igualdade por —1 e tomando t = -t temos -x(t) = 
x(-t) e y(-t) = y(t) =-y, mas x(t) = 0, logo x(-t) = O. Então x(-t) =0 e y(-t) = —y.. 


Figura 3.4 — Campo simétrico com relação ao eixo y que não tem órbitas periódicas. 


Fonte: Elaborado pelo autor. 


Vejamos agora alguns resultados de não-existência de órbitas periódicas. Seja F e 
C'(A, R) definido no aberto A c R2. 

Proposição 3.3.4 (Método de Poincaré das Curvas Tangenciais) Considere a 
família de curvas F (x, y) = C, com CE R. Se a função 
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of 


pp, 
[oh] ii 


OF 
“ay 
tem sinal constante em A e o conjunto de pontos de A nos quais 
oF oF 
P— — =0 
dr o Oy 


não contém uma órbita completa de (3.1), então o sistema (3.1) não possui uma 
órbita periódica que esteja inteiramente contida em À. 
Demonstração: Suponha que 


por QE o 2.4) 
dx Oy 
em A, 


Digamos que y é a imagem da aplicação t » (x(t), y()) de período T > 0. Se 
integramos F ao longo de y temos 


Tor ua 
fr=[ ECO NTO MORRE CIO O) A 
Ed 
E / SP (a(o). u(t))dt = Pl), (7) = Flo(0), (0) 


Como y é órbita periódica de período T, temos então que o lado direito da igualdade 


fr 


0= [= [o Ple(onS(r( U() + OE a) (ot) 


é nulo, logo, 


Mas daí, 


De (3.4) temos que o integrando da integral à direita é >0 em [0, T]. Logo, PoF/ox + 
QoF/oy = O sobre y, contradizendo a hipótese. 


Teorema 3.3.5 (Bendixson) Suponha que a função divergente OP/Ox + OQV/Oy do sis- 
tema (3.1) tem sinal constante em uma região simplesmente conexa R, e não é identica- 
mente nula em qualquer sub-região de R. Então o sistema (3.1) não tem órbita periódica 
que se encontra totalmente contida em R. 

Demonstração: Suponhamos que o sistema (3.1) possua uma órbita periódica y c 
R. Como R é simplesmente conexa, então y e seu interior estão totalmente contidos em R. 
Agora pela fórmula de Green temos 
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frio — Qda = /f (o + 3) dedy. (3.5) 


Mas Pdy — Qdx = O ao longo de toda y, dado que y é uma órbita. Assim temos que 
o lado esquerdo de (3.5) é zero, mas o integrando do lado direito tem sinal constante e 
não é identicamente nulo em S, logo a integral dupla é diferente de zero, o que gera uma 
contradição. 

Observação 3.3.6 Com o uso da fórmula de Green, segue que, em qualquer região 
limi- tada por órbitas periódicas, a integral da divergência, conforme definido em (3.5) é 
igual a zero. 

Teorema 3.3.7 (Dulac) Se existe uma função real B, de classe C!, em uma região 
sim- plesmente conexa R tal que BP )/Ox + YBQ)/0y de sinal constante e não nula em 
qualquer sub-região de R, então o sistema (3.1) não possui órbita periódica completa- 
mente contida em R. 

Demonstração: A prova é análoga a do Teorema 3.3.5, utilizando BPe QP no lugar 
de Pe Q, respectivamente. 

Observação 3.3.8 Chamamos uma função B como do Teorema 3.3.7 de função de 
Dulac B(x, y), e o método de prova de não existência de órbita periódica, dado no Teorema 
3.3.7, é chamado de método da função de Dulac. 

Os Teoremas 3.3.5 e 3.3.7 podem ser generalizados da seguinte forma. 

Teorema 3.3.9 (Dulac Generalizado) Se mudarmos a região R do Teorema 3.3.5 e 
no Teorema 8.3.7 por n-multiplamente conexa (isto é, R possui uma ou várias curvas exteri- 
ores limitadas e n — 1 curvas interiores limitadas), então o sistema (3.1) tem no máximo 
n— 1 órbitas periódicas que estão inteiramente contidas em R. 

Demonstração: A partir da prova do Teorema 3.3.5, sabemos que, se houver uma 
trajetória fechada y do sistema (3.1) em R, então y deve conter pelo menos uma curva 
limitada interna C de R em seu interior. Da mesma forma, também sabemos que se o 
interior de y também contém outras trajetórias fechadas y,, ..., y, então a região no interior 
de y, mas no exterior de todas as trajetórias y,, ..., Y, também contém pelo menos uma 
curva limitada interna C, por causa da Observação 3.3.8. Dizemos que C corresponde a y. 
Podemos ver que para diferentes y as curvas correspondentes C também são diferentes. 
Portanto, se o número de trajetórias fechadas em R for superior a n — 1, o número de 
conexos de R deve ser maior que n. Assim o teorema é provado. 

Daremos agora um exemplo da aplicabilidade do Teorema de Dulac. O exemplo 
pode ser encontrado em X. HUANG, Y. WANG and A. CHENG [15]. 

Exemplo 3.3.10 (Um sistema presa-predador cúbico) Seja o sistema 
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— = e(m + apt — aga”) — hay, 


dA l+r—y). 


Sea +a,>a,e k1,o sistema (3.6) não tem ciclo limite. 
Demonstração: Consideramos a função de Dulac B(x, y) = x?y!' e seja 


P(x,y=x(a + ax -— ax?) — kxy, x, )=y-1+x—)). 


Então, temos 


"9BP | 9BQ 


dvBR BQ)= "+53 
E y 


-a, xy! -ay!+(k- 1)x?. 


Se kz1, div(BP BQ) <0, e em qualquer sub-região do domínio, div(BP BQ) = 0. 
Então pelo Teorema 3.3.7, o sistema (3.6) não possui ciclos limites. 
Proposição 3.3.11 Se oP/ox + OQ/0y = 0, isto é, 


Pdy — Qdx = O (3.7) 


então o sistema (3.1) não tem ciclo indefinido, nem ciclo limite externo e nem ciclo 
limite interno em A. 

Demonstração: Seja a integral geral de (3.7) dada por &(x, y) = C que representa 
a equação de uma família de trajetórias do sistema (3.1). Ao longo de cada trajetória, & 
tem valor constante, e ao longo de trajetórias diferentes, mas adjacentes, & não pode 
ter o mesmo valor. Se houver um ciclo limite y de (3.1), então cada trajetória em uma 
vizinhança suficientemente pequena de y tem y como w(a)-limite. Daí pela continuidade 
de d(x, y) sabemos que nessa vizinhança D(x, y) = C, o que é impossível. Agora, se y 
é ciclo limite indefinido, então em trajetórias não fechadas que cobrem uma região anular 
aberta próxima de y na qual &(x, y) deve ser identicamente igual a uma constante. O que 
não acontece devido ao teorema de Green para órbitas constantes e por um lado só ter 
órbitas periódicas. 

Introduziremos a seguir resultados que permitem localizar órbitas periódicas de 
campos planares. 

Uma função Re C'(A, R) tal que 


papas , Yz,y)€A, (3.8) 
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é o fator integrante do campo de vetores X. A função R é um fator integrante de X 
em Ase, e somente se, R é a solução da equação diferencial parcial 


OR OR oP 0Q 
= 3.9 
o ay (+55) R (3.9) 


ema. 
A função V: A > R de classe C', é um fator integrante inverso do campo de vetores 
X, se Vverifica a equação diferencial parcial 
p? | a sh = (5 n 3) v 
[oh Oy dx Oy (3.10) 
em U. Notemos que V satisfaz (3.10) em A se, e somente se, R= 1/V satisfaz (3.9) 
em M(x )peA:Ví(x y=0 


Teorema 3.3.12 Seja (P Q) campo de vetores C! em um subconjunto aberto U de 
R2, esejaV=V(x, y) uma solução de (3.10). Se y é um ciclo limite do campo de vetores 
(P Q) no domínio simplesmente conexo de definição de V, então y está contido em > =((x, 
Wjeu:Ví(x,y)=0). 

Demonstração: A existência do fator integrante inverso V, definido em uma região 
simplesmente conexa U implica que o campo de vetores (P/V, Q/V) é hamiltoniano em U 
>. Uma vez que o fluxo de um campo vetorial hamiltoniano preserva a área, e em uma 
vizinhança de um ciclo limite, um fluxo não preserva a área, segue o teorema 


No Capítulo 4, veremos uma aplicação do Teorema 3.3.12. 
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O TEOREMA DO "AVERAGING" CLÁSSICO 


O método do 'Averaging' é uma ferramenta para o estudo de sistemas não lineares 
com coeficientes periódicos. O método consiste essencialmente em associar ao sistema de 
EDO's em estudo, um novo sistema chamado de sistema promediado, de modo que cada 
singularidade do sistema promediado esteja associada a uma órbita periódica do sistema 
original. 

Neste capítulo trataremos apenas do caso do 'Averaging' periódico de primeira or- 
dem. A primeira formalização dessa teoria foi feita em 1928 por Pierre Joseph Louis Fatou. 
Este capítulo da dissertação está baseado em S. SCHECTER and F. SINGER [14], J.A. 
SANDERS, F. VERHULST and J. MURDOCK [04] e F. VERHULST [02]. 

Este capítulo está organizado do seguinte modo. Na Seção 4.1 é feita a apresen- 
tação do método do 'Averaging' clássico de primeira ordem. Os resultados preliminares a 
sua demonstração estão na Seção 4.2. A demonstração do método é feita na Seção 4.3, 
enquanto que na Seção 4.4 aplicamos o método para estudar a existência e estabilidade 
de um ciclo limite para a equação de Van der Pol. 


4.1 O MÉTODO DO “AVERAGING' PARA CAMPOS DE VETORES SUAVES 


Consideremos o sistema diferencial 


Ú =eF(ta) + R(t,2,e), 


x(0) = xo. (4.1) 


com xe Dc Rº, D um conjunto aberto, limitado e conexo, e t > 0. Além disso, 
assumimos que F: [0,+0)x D> Rre R:[0,+0) x Dx (-e, e) — R" são funções suficien- 
temente regulares, sendo e, > 0. Suponha ainda que Fe R'são T -periódicas na primeira 
variável, isto é, existe T>0Otalque F(t+ Tx) =F(t x)e A(t+ T, x, e) = R(t x, €), para todo 
te [0, +00), xe De e E (-€,, €,). No momento apropriado estabeleceremos a regularidade 
de Fe R. 

Simplificando esse sistema através da eliminação do termo de ordem e? para que o 
campo se torne autônomo, calcula-se a média com respeito a t do sistema. Então, a ideia 
do método é considerar a velocidade média do fluxo, pois assim, podemos definir um novo 
sistema que é autônomo e então comparar as respectivas soluções. 

O sistema promediada associado ao sistema (4.1) é definido por 


p, — efO(» 
“0) A ' me) 
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onde a função fº: D > R" é dada por 
i 7 
fz) = DI) F (t, z)at. (4.3) 
TJ 


Chamamos fº de função promediado associada a (4.2). 

O principal objetivo é obter informações sobre existência e estabilidade de órbitas 
periódicas do sistema (4.1) através das singularidades do sistema (4.2). Mais precisa- 
mente, enunciamos a seguir o principal teorema deste capítulo. 


Teorema 4.1.1 ((Averaging' de primeira ordem clássico) Consideremos o sistema: 


x(t) = eF(t x) + ER(t x, €), (4.4) 
comt>0,xe D, Dc R" um conjunto aberto, F:[0,+0)x D> RreR:[0,+0)x D 
x (-€,, 6) > R”, sendo e, >0. Suponha que: 


1. as funções F, R, OF/dx, PF/oxº e OR/Ox são contínuas e limitadas por uma cons- 
tante M>0; 


2. as funções Fe R são T-periódicas na primeira variável. 


Então são válidas as seguintes afirmações: 


a. sepé um zero da função promediada (4.3) tal que 


of? 


então existe uma solução T-periódica t » x(t, e) do sistema (4.4) tal que xX(0, > p 


=p 0, (4.5) 


quando € > 0; 


b. se todos os autovalores de 


(5) 
dz ) (4.6) 


possuem parte real negativa então, para lel >0 suficientemente pequeno, a solução 
periódica t » x(t, e) correspondente do sistema (4.4) é assintoticamente estável, e se um 
dos autovalores de (4.6) possui parte real positiva então t» x(t, e) é instável. 

Aqui o símbolo dF/ox denota a matriz n x n formada pelas n últimas colunas da matriz 
jacobiana de F, análogo para dR/0x. O símbolo &F/dx denota a matriz hessiana de F. 

A ideia da prova do Teorema 4.1.1 consiste em introduzir uma mudança de coor- 
denadas próxima da identidade que é usada para simplificar a equação (4.1). A hipótese 
(4.5) é utilizada para aplicar o Teorema da Função Implícita, o qual nos remete a uma 
vizinhança onde existe uma solução periódica de (4.1). 
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Antes da demonstração do Teorema 4.1.1, vejamos alguns resultados de caráter 
técnico que nos auxiliarão nessa tarefa. 


4.2 RESULTADOS ÚTEIS 


Os próximos resultados serão utilizados na demonstração do Teorema 4.1.1. Aque- 
les que consideramos mais conhecidos terão suas demonstrações omitidas. 


Teorema 4.2.1 Seja fcontínua no aberto A deR x Ex 4, com EC R, Aaberto de Rº 
e D,f contínua em A. Então, para À fixo, a solução q = q(t, ty, Xy À) de 


T=efttxi ME, 


0 


é única e admite derivada parcial D,f com relação a x,. Mais ainda, a aplicação (t, t, 


Xy À) > D,op(t to, X, À) é contínua no seu domínio 


D=(tt,X, Ni (tt,Xx, NEM 0 (t, x, <t<o (t, x, ND 


xX)=D,o(tt,X, A) e,= de (é to Xy À), 
xy. 
para todo 1 < k < dimE, é solução de 
x =x, x(t) = e, 
onde 
HD = Ht to x, A) =Df(t (tt X À), À). 


1 tor 


A demonstração pode ser encontrada em J. SOTOMAYOR [09]. 


Teorema 4.2.2 Consideremos as funções F,R e fº dadas nos sistemas (4.1) e (4.2). 
Suponhamos que F é lipschitziana na segunda variável e R é contínua. Então existe uma 
constante C >0 tal que 


É 


lIx(t e) — Z(t elll< Ce para0O<se<se e0O<t< 
E 


onde x(t, €) e Z(t, €) são soluções de (4.1) e (4.2) respectivamente. 
A prova do Teorema 4.1.1 que apresentaremos aqui, utiliza uma mudança de coor- 


denadas próxima da identidade que transforma o sistema (4.4) em um sistema da forma 


Z=efAD)+e st ze), 
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sendo S uma função apropriada. Para isso precisamos do resultado que segue, que 
pode ser encontrado em J.A. SANDERS, F. VERHULST and J. MURDOCK [04]. 

Teorema 4.2.3 (Perturbação da Identidade) Considere um aberto DC R"e e >0. 
Sejay: Rx Dx[0, e] > R" uma aplicação de classe C' e lipschitziana na segunda variável, 
isto é, existe M, >0 tal que 


It, Y, 6) — ut Y» Olls My, -YlLteR,y, y,€D,celo, e]. 


Então existe e, >0 tal que, para cadate Re ee O, e] fixados, a aplicação 
Ultye=y+entye,yeD, 
tem uma inversa da forma 
w(t x e)J=x+evtt, x,e), 
definida em um conjunto U (t, D, e) c Rº, sendo v uma aplicação de classe C!. 
Demonstração: Escolha e, talque O <e, <min(e, 1/M >. Começamos demonstrando 
a injetividade de U para t E Re e, € [0, e] fixados. Dados y,, y, E Dcom y, * y,, temos 
HU (ti, Y, E) — Ult, Yo el =Iy, — y, — Mt, Yo E) — Mto Y, ED 
2 ly, — yjll — ellu(t, Y,» E) — Mt, Y, El 
= ly, — yl — eMly, — yj 
=(1- SoM y, = Yall. 
Como 0 < «,M, <1 segue a injetividade. Logo, a função U tem inversa única definida 
em um conjunto U(t, D e) c Rr. 
Utilizaremos o Teorema da Função Implícita para demonstrar a regularidade e a 
forma da função inversa. Considere a aplicação 


V(tyexy=x-y-etye),teR, yeD,ecelO, e] xeR”. 


Para cada y, € D, tome x, = y, + €,V(t, Yy €,), temos então que V'(t, Y, 6 X) =0. 


Como 
oV 


Em (b, Yy O, X) =—1, 


segue do Teorema da Função Implícita que existem vizinhanças V, e |, x J,x W,, de 
Ye (t, O, x) respectivamente, e uma aplicação h: | x J,x W, > V, de classe C', tais que 


x-h(te,x)-eut ht ex), J=0, (Le M)MElxJxW, 
obtendo assim 


h(texM)=x- et he x, e), (beM)EelxdxW. 
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É claro que podemos supor que Sc lo, e]. 

Como já demonstramos que a função U tem uma inversa global, esta coincide com 
cada uma das inversas locais obtidas em seus respectivos domínios. 

Os próximos resultados serão úteis na demonstração do item (b) do Teorema 4.1.1. 

Lema 4.2.4 (Desigualdade de Gronwall) Sejam q, y :[t, t,+ a] > R funções conti- 


nuas e não-negativas definidas no intervalo [t + al, sendot eRea>0 constantes. 


o bo 


Suponha que 
)<o, y W(s)p(s)ds+ 0, telt, t+ al, (4.7) 
onde 6,>0 para i=1,2. Então 


o(b < ethos felt, t+ al. 

A demonstração pode ser vista em J. K. HALE [06]. 

Estabelecemos a seguir o decaimento exponencial das soluções de sistemas lIi- 
neares com coeficientes constantes, no caso em que todos os autovalores da matriz do 
sistema têm parte real negativa. 

Proposição 4.2.5 Considere o sistema diferencial linear 

T=aAsx, (4.8) 
onde A é uma matriz real constante n x n com autovalores À,, ...A,. Então valem as 
seguintes afirmações: 
1. (a)seReA,<0 parak=1,..., n, então para cada solução x de (4.8) tal que x(t,) 


=X, comtpe Rex,e€ Rº arbitrários, existem duas constantes positivas C eu 
satisfazendo 


IIx(BI < Clbg let, t>0. 


2. (b)seReA,<0 parak=1,..., n, então toda solução x é limitada para tz 0. Mais 
precisamente existe C >0 tal que 


IIx(BIl < Clix Il, tt; 


2 lo 


3. (c) se existe um autovalor À, com Re À, >0, então para cada vizinhança V da 
origem, existe uma solução x do sistema (4.8) tal que xX(0) E Ve 


lim Ilx(BIl = 00. 
t>+oo 
Seguindo a noção de estabilidade no sentido de Liapunov, na Proposição 4.2.5 
temos que soluções como as dos itens (a) e (b) são casos especiais de soluções assin- 
toticamente estáveis e estáveis, respectivamente. Uma solução como na afirmação (c) é 


chamada de instável. 
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O próximo resultado também é conhecido como Teorema de Poincaré-Liapunov, e 
estende a Proposição 4.2.5 localmente para sistema não lineares. 

Teorema 4.2.6 Sejam D um aberto de R" contendo a origem, | um intervalo aberto 
deRete |, x, € D arbitrários. Considere o sistema diferencial 

x = Ax+ B(bx+f(t x), x(t) =X, (4.9) 

onde te R, A é uma matriz real constante n x n com todos os seus autovalores 
com parte real negativa, B é uma matriz real n x n cujas entradas são contínuas em | tal 
que lim,.. IB(BIl = O. Suponha que a função f seja contínua e lipschitziana em x numa 
vizinhança de x= 0. Se 


lim Ft, m) 


= O uniforme em t, 
lizilo |] 


então existem constantes positivas C, d ey tais que Ilx ll sô implica 
Ix(BIl < Clix leio para t > t,. 
Apresentaremos a seguir o Teorema de Floquet e algumas de suas consequências. 
Maiores detalhes podem ser encontrados em J. K. HALE [06]. Consideremos agora o 
sistema linear diferencial da forma 


x = aA(tx, (4.10) 


onde A é uma matriz real n x n, com entradas contínuas em Re T -periódicas. 

Observe que mesmo A sendo periódica, as soluções do sistema (4.10) podem não 
manter essa propriedade, como mostra o exemplo %: = (1 + sen 2x, cujas soluções são da 
forma x(t) = celcost, 

O próximo resultado, geralmente chamado Teorema de Floquet, diz que a matriz 
fundamental do sistema (4.10) pode ser escrita como um produto de uma matriz T-periódica 
e uma matriz não periódica. 

Teorema 4.2.7 Considere o sistema linear diferencial (4.10) sendo A uma matriz n 
x n, com entradas contínuas e T-periódicas. Então, cada matriz fundamental & do sistema 
(4.10) pode ser escrita como um produto de duas matrizes nx n 

o()=P(De, te R, 

onde P é T-periódica e B é uma matriz constante. 

O próximo resultado permite obter informações do sistema (4.10) a partir de um 
sistema de coeficientes constantes. 

Corolário 4.2.8 A mudança de coordenadas x » P (t)x transforma o sistema (4.10) 
no sistema de coeficientes constantes 

T =Bx. 
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Definição 4.2.9 (Matriz de monodromia) Uma matriz de monodromia do sistema 
(4.10) é qualquer matriz não singular C tal que P(t+ T) = AC, te R, sendo & uma matriz 
fundamental de (4.10). 

Definição 4.2.10 (Multiplicadores característicos) Os multiplicadores característicos 
são os autovalores de uma matriz de monodromia C. 

Observação 4.2.11 Se d e 'P são matrizes fundamentais de (4.10) e D(t+ T) = O(f) 
C, para qualquer te R, então existe matriz não singular D tal que 'P(?) = P(DD. Logo P(t + 
T)=O(t+ T)D= O()CD = 'W()D'CD e assim D'CD é matriz monodromia relativamente 
ap(d. 

Assim as matrizes de monodromia relativas a matrizes fundamentais distintas são 
similares, possuindo então os mesmos autovalores. Concluímos que os multiplicadores 
característicos são unicamente determinados. 

Observação 4.2.12 (Expoente Característico) Um expoente característico é 
qualquer número complexo À tal que eN é multiplicador característico. 

Verifica-se que sempre podemos escolher os expoentes característicos como os 
autovalores de uma matriz B satisfazendo C = e?º. 

Apresentamos a seguir uma versão do Teorema 4.2.6 para sistemas periódicos. O 
resultado pode ser encontrado em J. K. HALE [06]. 

Proposição 4.2.13 Considere o sistema diferencial 


x =aA(tx+ ft x), (4.11) 


em que A é uma matriz real nxn com entradas contínuas e T-periódicas, f: RxR” 
> Rn contínua na primeira variável e numa vizinhança da origem na segunda variável. 


Suponha que 


lim Ft) 


= O uniforme em t. 
Iell5o |[x|| 


Se a parte real de todos os expoentes característicos do sistema 
y = A(Dy (4.12) 


são negativos, então a solução x = O do sistema (4.11) é assintoticamente estável. 


4.3 DEMONSTRAÇÃO DO TEOREMA DO 'AVERAGING' 


Finalmente, passamos à prova do Teorema 4.1.1. 
Demonstração: (Afirmação (a)). 
Introduzimos a função T -periódica: 
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u(t x) = / (F(s, x) -fM)ds, tz0,xe D, 


Afirmamos que 
Iu(t Ml<2MT t>0,xe D, (4.13) 


onde Me T são dados no enunciado do Teorema 4.1.1. 
De fato, para0<t<sTexe Dtemos: 


inteali= | [riso -1taas 


/ F(s,x)ds — / fº(a)ds 


/ F(s,x)ds Ih f(a)ds 


< [ Plso)tas+ [Il (o)lds 


to t ' | 
<[ Mas+ | If (a)lds. 


< — 


Contudo, 


T T 
< a) F(s,a)ljds < z) Mds = =M. 


1$º(a)|| = Ef risos 
Assim, 


t t 
Hu(t, x) <2[ Mds = 2m | ds, 
0 0 


logo 
Iu(t, Mll<s2MT, O<t<T (4.14) 
Se t>T,temos que existe me Ntalque0 <t-mT=<T. Logo, 
t 


m iT 
ue=5 [o (Pl) 1o)as+ [ (P(s,0) = Pos 


j=i (G—-DT mT 


Como Fe fº são T periódicas, para j= 1,2, ..., m, segue que: 
Er t-mT 
u(t X) = m[ (F(s, x) - f(x))ds À (F(s, x) — fO9))ds. 
0 0 


Da definição de fº resulta que a primeira integral do lado direito da igualdade acima 


é zero. Logo 
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u(t %) = / Cedo (4.15) 
(0) 


provando a afirmação, por ser reduzida ao caso em que O <t< T. De (4.14) segue 
então (4.13). 
Afirmamos agora que 


(t, 2) 


[ao ti )|s2MT; tER, XE D. (4.16) 


De fato, para cada j= 1,2,..., nte Re xe Dtemos 


| Deli 0] E EMP — fº(x))ds 
= Ed z) Seta as) 
E 0 Oz; RE 


OF á 
dx 8,5) Ox (2) 


ds 


<[ 
0 


t 0 
a] E) ds + [E (a))|l d 
0 Lj 0 
logo 
9) t 9) 0 
[De(s <MT+ | a ds. (4.17) 
dz; 0 Ox; 
Mas, 
of “jo (1 r o 
à] no Fls 
rf. : of ,a)ds 
si. do (82) ds<M. 
Portanto, de (4.17) segue que 
Õ 
|i(s.a) <MT+M. 
OE; 


Como podemos supor que T > 1, segue (4.16). 

Diminuindo D se necessário, podemos supor que p € De que D é um aberto conexo 
com fecho compacto. De (4.16) segue que podemos aplicar o Teorema do Valor Médio e 
concluir que uy é uma função lipschitziana na segunda variável. Assim, pelo Teorema 4.2.3, 
estabelecemos a mudança de variável: 
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x(D) = Z(t) + ext, Z(t)), (4.18) 


no qual € E [0, € ] e e, é dado no Teorema 4.2.3. Reescreveremos o sistema (4.4) 
na variável z. Para simplificar a notação omitiremos a dependência de x e z da variável t. 
Assim, diferenciando (4.18) em relação a te utilizando a equação (4.4) obtemos a equação: 


2+e-(t)re—,2)2=t=<F(t,s)+R(,2,0), 


[h — a, 2) 2=eF(t,2) + é R(t,2,0 — co, z). (4.19) 


Da definição de y segue que 


de 2=F(t 2) — f(2), 


logo de (4.19) temos 


ou E 
+ Ema lt» z)| Z=eF(t x) + rt x, 6) - elF (t 2) - (2). (4.20) 
z 
Seja 
S=eF(tz+edt 2))-eF(t 27 + eR(t z+ et 2),€), (4.21) 


então podemos reescrever (4.20) como: 


[ + Rm 2] Z=(7+5S. (4.22) 
oz 


Como &F/o0xº é limitada e D é conexo com fecho compacto, segue que dy/Oz é 
lipschitziana na segunda variável. Assim, diminuindo e, se necessário, pelo Teorema 4.2.3 
podemos inverter a matriz do lado esquerdo de (4.22). Da fórmula 


(L+AJ!=| —-A+A-.., 
válida para matrizes tais que IIAIl < 1, obtemos então a seguinte expressão: 


[1 + co, gi t=doo A o(e), (4.23) 


paratz0ezeD. 
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Afirmamos que, para S definida em (4.21) temos S = O(e?). De fato, como oF/ox é 
limitada e D é conexo, do Teorema do Valor Médio segue que F é uma função lipschitziana 
na segunda variável. Denotamos a constante de Lipschitz por L, >0. Temos que 

ISll=IleF(tz+eMt z2)-eF(tz)+e AR(t z+ eut 2), elll 
seL lui 21 + e2lA(t z+eu(t 2), el. 


Da estimativa (4.13) e da hipótese 1 do Teorema 4.1.1, 
IS <eLp att 2) + ARE z+ ex(t, 2), e)ll 
s el, 2MT + eM 
< ek, 
sendo K definido por K= L, 2MT + M. Assim, 


S = O(e?). (4.24) 
Portanto, de (4.22), (4.23) e de (4.24) temos que Z' satisfaz uma equação da forma: 
s ou 
z=ef(7)+S-€ dz 22) + Ole). (4.25) 
Utilizando a Fórmula de Taylor, temos que 


F(tz+eMt 2)-Fl(tz)=e EM zZ) + O(e?) (4.26) 


= Put, 2) + O(e?. 


Definamos agora S (t, z, €) por: 


S(t,z,e) = >—(t, z)u(t, 2) — o (2) + R(t, 2,0) +0(6), 


onde A(t z + eu(t 2), €) = R(t z, e) + O(e), que é obtido através da Fórmula de 
Taylor de primeira ordem. Combinando (4.26), (4.21) e A(t, z+ ext, 2), e) = A(t, z, e) + Ofe), 
podemos escrever (4.27) como: 


s=A7+esS(t z 6). (4.27) 
Observemos que O(eº) = e“O(e), e portanto, podemos incorporar tal símbolo no 
termo est z, 6). 


Pelas condições de regularidade nas hipóteses, temos que S(t, z, €) é de classe C! 
com respeito a z, e que é T -periódica pois //, F, Re o erro O(€*) são T -periódicos. 
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Mostraremos a seguir que, pela escolha de y, uma solução T -periódica z de (4.27) 
corresponde a uma solução T -periódica x de (4.4). De fato, suponhamos que Z(t) = Z(t, + 
T), então 


x(t) — ext, At) =x(b+ T)- eti+ T t+ T)). 
Mas 


to+T 
guto + P,z(to + T)) E [F(s, z(to + T)) — fº(z(to + T))lds 


[=] 


— $(z(to)lds 


— 


to “tp 


to)) — fº(z(to))lds 


+ 
+ 


(s, 2(to)) — fº(z(to))lds 


Ta 


+ 


[F(s, 2(to)) — $(z(to))lds 


“rt, Ato)) — fº(c(to))lds 
to, Ato). 


Es 


Assim x(t) = x(t, + T), e portanto x é T -periódica. 
Escrevendo (4.27) na sua forma integral temos: 


z(t) + f(«s))ds + e “36, z(s), eds. 


Como uma solução z da equação acima é T -periódica se, e somente se, Z(0) = Z(T), 
temos que uma solução 7 -periódica z do sistema (4.27) é solução de: 


T 
frias) (spas +e [ S(s, Z(s), eJds = 0. (4.28) 
0 


Considere a aplicação w = w(t, a, e) tal que ww é a solução do sistema (4.27) com 
w(o, a, e)= a, onde ae De e E [0,€,) são quaisquer. Como as soluções de (4.27) dependem 
diferenciavelmente dos parâmetros (veja o Teorema 4.2.1), segue que w é de classe C!. 


Observe que se e = O então de (4.25) resulta que w é a função constante igual a a, logo 


O rapel; teR,a ED. 
A 


Defina 


T Es 
h(a, e) = / Fw(s, a, e))ds + ef S(s, w(s, a, e), ejds, 
b 0 
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então = 
E ) f(pjds 
0 


Por hipótese fº(p) = 0, logo h(p, 0) = O. Além disso, 


é (a,0) = e 1 (w(s,a,0))ds 
= 2ru(s, a,0))ds 
= = él a 0) Se(s, a, O)ds 
= 2 f(w(s,a,0))Inds. 
Logo, 
Son, 0) = TES) 

Da hipótese 2 do Teorema 4.1.1 segue que dh/da(p, 0) é inversível. 

Pelo Teorema da Função Implícita, existe um aberto Vc D, tal que (p,0) E Vx [0, 
€,) e, para cada E E TO, €,), existe um único a E V'tal que a solução do sistema (4.27), com 
Z(0) = a, satisfaz (4.28). 

Ainda pelo Teorema da Função Implícita, temos que Z é uma função diferenciável 
com relação ao parâmetro e. Em particular, da sua continuidade, podemos inferir que Z(e) 
> psee>0. 

Demonstração da (Afirmação (b)) Começamos linearizando a equação (4.4) em 
uma vizinhança da solução periódica t » x(t e). Depois, fazendo uma translação x = z 
+ x(t, €), expandindo com respeito a z, omitindo os termos não lineares e renomeando 
as variáveis dependentes novamente por x, obtemos a equação diferencial linear com 
coeficientes T -periódicos 
T =eA(t e)x, (4.29) 


onde 
A(t, €) = SE, x) + eR(t x, €)] 


xexeO 


Definimos a matriz T -periódica 


B(o) = So(t,0) 


e pelo item (a) temos que lim,  A(t e) = B(?). Definamos também as matrizes 


1 Y g o T = ] 
=) B(t)dt e c(m=[ [B(s) — BºJds. 


Note que Bº é a matriz promediada do sistema linearizado. A matriz C(t) é T- periódica 
e a sua média é zero. Tomando a transformação definida por y = (| — e())x obtemos 
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y =eC (x+(1- Di 

— elx + eBº x(l- eC(D)eA(t e)x (4.30) 
[eBº + e(A(t e) — B(9)) - eC(DA(L e(!- ec(d)!y 

eB'y + e(A(t e) — B(D)y + eS(t e)y. 


A função S(t, €) é T-periódica e limitada. Notemos que A(t, e) — B(t) > O quando 
€ > 0, e também que o expoente característico do sistema diferencial (4.30) depende 
continuamente do parâmetro e. Portanto, para e suficientemente pequeno, o sinal da 
parte real dos expoentes característicos é igual ao sinal da parte real dos autovalores da 
matriz Bº. Obtemos a mesma conclusão, usando a transformação y = (| — eC(t))x, para os 
expoentes característicos do sistema diferencial (4.29). 

Aplicando agora a Proposição 4.2.13, obtemos a estabilidade da solução periódica 
no caso em que todos os autovalores possuem parte real negativa. Se um autovalor pos- 
suir parte real positiva, a transformação de Floquet e a aplicação da Proposição 4.2.7 
fornecem a instabilidade da solução periódica. 


4.4 EQUAÇÃO DE VAN DER POL 


Apresentaremos uma interessante aplicação do Teorema 4.1.1 a equação diferen- 
cial de Van der Pol 
T+x=e1-3T, 
a qual pode ser escrita como campo de vetores suave como 
T=y, 
y=-x+e(1 — X)y. (4.31) 


Em coordenadas polares (r, 0), onde x=rcos 0, y=r senô, esse sistema torna-se 


r=er1- cos? 6) sen?6, 
6=-1+€6(1 - P cos? 8) sen?o. (4.32) 


Agora tomando 6 como nova variável independente em (4.32) e expandindo em 
série de Taylor em e obtemos 


— =- er(1 - P cos? 60)zsen?0 +0(e”). 


Notemos que o campo de vetores acima está na forma padrão para aplicar o Teorema 
4.1.1. Tomando x=r, t=0, T=2mne F(t x) =-r(1 — ? cos? 6) sen?0, temos que 
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e 1 
F!(r) = >|, r(1 — r2 cos? 8) senZ6 dg — dd —4). 


A única raiz positiva de F '(1) é r= 2. Como dF '/dr(2) = 1, pela afirmação (a) do 
Teorema 4.1.1, segue que o campo (4.31) tem, para lel * O suficientemente pequeno, um 
ciclo limite bifurcando da órbita periódica de raio 2 do sistema não perturbado (4.31), com 
e = 0. Além disso, como dF '/dr(2) = 1 >0, pela afirmação (b) do Teorema 4.1.1, este ciclo 
limite é instável. 
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REALIZAÇÃO DE UMA CONFIGURAÇÃO DE CICLOS LIMITES 


Neste capítulo mostraremos uma espécie de recíproca do XVI problema de Hilbert, 
isto é, dado um conjunto finito qualquer de curvas de Jordan no plano, duas a duas disjuntas, 
é possível exibir explicitamente um sistema de EDO's polinomial, que tem exatamente tais 
curvas como ciclos limites. Mais do que isso, o resultado estabelece uma relação entre o 
número de curvas de Jordan dadas e o grau do polinômio que define o sistema de EDO's. 

A única ressalva feita sobre o estudo do teorema central deste capítulo é de que não 
temos controle sobre a estabilidade dos ciclos limites construídos, pois podem ser estáveis, 
semi-estáveis ou instáveis. 


5.1 O XVI PROBLEMA DE HILBERT 


O XVI problema de Hilbert ainda é um dos poucos problemas em aberto da famosa 
lista de problemas formulada pelo matemático alemão David Hilbert durante o Congresso 
Internacional de Matemática de Paris em 1900, lista qual influenciou, e influencia ainda, 
muitos matemáticos que buscam por seu nome na história dos números. Ele trata dos 
ciclos limites em sistemas diferenciais polinomiais e pode ser dividido em duas partes. A 
segunda parte do problema tem o seguinte enunciado: qual o número máximo de ciclos 
limites de um sistema diferencial polinomial, e o que podemos dizer sobre as posições 
relativas dos seus ciclos limites? 

J. ÉCALLE [05] e Y. I"YASHENKO [16] mostraram que o número de ciclos limites 
de um sistema polinomial planar é sempre finito, mas não existem resultados que mostrem 
uma cota uniforme para esse número. Em contrapartida, existem exemplos mostrando que 
se o campo não for polinomial, o número de ciclos limites pode ser infinito. 


5.2 O TEOREMA PRINCIPAL 


O enunciado preciso do resultado principal deste trabalho requer algumas defini- 
ções preliminares. 
Definição 5.2.1 


a. Uma configuração de curvas fechadas é um conjunto finito C=(C,, .... CJ de 
curvas fechadas, simples e, duas a duas disjuntas do plano. 


b. Dada uma configuração de curvas fechadas C =(C,, .... C), a curva C, é dita 
primária se não existe curva C, da configuração C contida em IntC, 


c. Duas configurações de curvas fechadas C=(C,,... Che C'=(C,,..., Ci) são 
topologicamente equivalentes se existir um homeomorfismo h : Rê > R2 tal que 
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h (U c) = U e 
1=1 i=1 


d. Um campo vetorial planar X realiza a configuração de curvas fechadas C como 
uma configuração de ciclos limites se o conjunto de todos os ciclos limites de X 
for topologica- mente equivalente a C. Neste caso, dizemos que a configuração 
C é realizável por X. 


e. Umciclo limite é dito algébrico se for uma componente de uma curva algébrica. 


Com a nomenclatura acima temos o seguinte resultado, devido a J. LLIBRE and G. 
RODRÍGUEZ [07]. 

Teorema 5.2.2 Sejam C=(C,, ..., C) uma configuração de curvas fechadas e ro 
número de curvas primárias de C. Então valem as afirmações: 


a. a configuração C é realizável por um campo vetorial planar polinomial, 


b. aconfiguração C é realizável como ciclo limite algébrico por um campo vetorial 
polinomial de grau menor ou igual a 2(n + 1) — 1. 

A afirmação (a) do Teorema 5.2.2 foi obtida pela primeira vez por S. SCHECTER 
and F. SINGER [14], mas nesse trabalho os autores não fornecem um campo vetorial poli- 
nomial explícito satisfazendo a dada configuração de ciclos limites, tampouco estabelecem 
a cota apresentada no item (b) do Teorema 5.2.2. Na demonstração do Teorema 5.2.2 
será construído explicitamente um campo vetorial polinomial de grau menor ou igual a 
2(n + 1) — 1 realizando C por ciclos limites algébricos. A ideia do modo como é feita a 
construção explícita do campo de vetores na demonstração do Teorema 5.2.2 provém de 
M. |. A!MUKHAMEDOV [12], que chama a construção de segundo problema inverso da 
teoria qualitativa das equações diferenciais, mas a referida demonstração tem algumas su- 
tilezas em contraponto com a forma apresentada em M. |. AL/MUKHAMEDOVY [12], o que 
permitem a construção explícita do campo de vetores polinomial. 


5.3 PRELIMINARES 


A prova do Teorema 5.2.2 requer alguns preliminares que apresentamos a seguir. 

Observação 5.3.1 Observemos que se C=(C,.., C)eC'=(C,,..., Ci) são 
configurações de curvas fechadas topologicamente equivalentes então n = m. 

Apresentaremos a seguir um exemplo de duas configurações de curvas fechadas 
que não são topologicamente equivalentes, apesar de serem formadas pelo mesmo nú- 
mero de curvas fechadas. 

Exemplo 5.3.2 As configurações de curvas fechadas a =(C,, C) e B=(C,, Cj) 
dadas na Figura 4.1 não são topologicamente equivalentes. 
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Figura 5.1 — Configurações a e P. 


B 


Fonte: Elaborado pelo autor. 


De fato, suponhamos que as configurações a e 8 são topologicamente equivalen- 
tes. Tome C, CiniC, p, EintCNC,U IntC,) e p, E IntC, [a -Figura 1.1], assim o segmento 
P1D2, intersecciona C,, pelo Teorema da Alfândega. Então, supondo que a e B são 


topologicamente equivalentes, temos que o segmento p/p), com p E IniC e p, E IniC,, 
deve também interseccionar apenas C,. Mas temos que, como por hipotese, Ce C, são 
disjuntas, e C, c ExtC,, então Disi(C,,C,) = €e>0 logo, o segmento PiDh, deve intersectar 
Ce C,,o que é um absurdo, pois o segmento não deve intersectar C, em . 

Nosso próximo objetivo é demonstrar que, para efeito de realização de uma dada 
configuração de curvas fechadas C, sempre podemos considerar os elementos de C como 
círculos. Para isso necessitamos do: 

Teorema 5.3.3 Dada uma configuração de curvas fechadas C, existe uma 
configuração de curvas fechadas C' =(C,, C;,..., CY topologicamente equivalente a C tal 
que C,é um círculo, para todo j= 1,2,..., n. 

A prova será feita usando uma versão do Teorema de Schonflies para grafos e 
apresentada no Apêndice. 

O teorema a seguir se deve a H. GIACOMINI, J. LLIBRE, and M. VIANO [03] e 
demonstra que os ciclos limites estão localizados nas curvas de nível de fatores integrantes 
inversos, quando este existe. 


Teorema 5.3.4 Seja X um campo de vetores C'! definido no aberto Uc Rº e seja V: 
U > R um fator integrante inverso de X. Se y é um ciclo limite de X, então y está contido 
emz=((xyjeuU:V(x, y)=0> 

Demonstração: Por hipótese, existe um fator integrante inverso V definido em U. 
Então, X/V é Hamiltoniano em UN. Então X/V e X são topologicamente equivalentes. 
Sabemos que o fluxo de um campo de vetores Hamiltoniano preserva área (o Teorema 
de Liouville garante que se um campo de classe C' tem divergência nula, então seu fluxo 
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preserva volume) e em uma vizinhança de um ciclo limite a área não é preservada, assim, 
temos que o ciclo limite y é UN, portanto, y c > 


5.4 DEMONSTRAÇÃO DO TEOREMA 5.2.2 


Seja C =(C, ..., C) a configuração de curvas fechadas dada no enunciado do 
Teorema 5.2.2. Do Teorema 5.3.3 podemos supor que cada curva C, j=1,..., n, é dada 
pela circunferência definida por 


fo y)=(X-x2+(y-y)2-1=0, (5.1) 


sendo (x, y) E RÊe 1, >0 são escolhidos de modo conveniente. Renomeando as 
curvas, podemos supor que Cé uma curva primária, para [= 1,2,...,r. 

Para cada curva primária C, j=1,..., r, definimos p= (X, y) e as funções com- 
plexas dadas por 


Pra ( y) = (x— x) + Ãy= y) 
Folk V) = (x— x) — ly — y). (5.2) 


Por conveniência, no momento apropriado faremos uma restrição no domínio das 
funções acima para R? — (1, onde Ó = Uj-i(-oo, x] x 0). 

Definamos agora o campo de vetores X abaixo e provaremos que satisfaz a condi- 
ção (b) do Teorema 5.2.2. Seja 


onde 
n+2r n+2r 
9fulx,y) 
Pt =-Sox (TT nim | He, (ee 
h=1 l=1Zk 
e 
n+2r n+2r 
O fula,y 
(e) = 55 M H TE aids. n) die (ER? 
k=1 [= 1k 
comA,=...=A,=1.Alémdisso, À, ,,,=1 +ieA,.,=1-iparaj=1,..,r 


Mostraremos a seguir que Pe Q são polinômios reais de grau <2(n + 1) - 1. Faremos 
a prova para P, e de forma análoga é mostrado que Q é polinômio real. Reescrevendo P 
da forma 
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n n+2r n+2r n+2r 
- (1 n)u- 5 1 1 
k=1 


[= lZk k=n4100 I=1%k 
temos 


(Li E (o) 


I=11Zk l=n+1 k=n+41 l=n+HIl£k 


(5.3) 


No primeiro somatório de (5.3), para k = 1, 2, ..., n, de (5.1) temos queTI;., 14h é 


Ra f 


um polinômio real de grau 2(n — 1), enquanto que mn n+1h 


é um polinômio real de grau 2r, 


pois 
zz =IH2 vzeEC. (5.4) 


Como 2y é um polinômio de grau 1, o primeiro somatório em (5.3) é um polinômio 
realdegrau=2(n-1)+2r+1=2(n+n-1. 

Denotamos o segundo somatório de (5.3) por S. Com o objetivo de analisarmos S, 
para cada k = 1,3,...., 2r— 1, definimos: 


n+2r—1 


EL fi 


I=n+1, 
l£n+k 
I£n+k+1 


que é um polinômio real de grau 2(r — 1) por (5.4). 


Então 
n+2r n+2r 
S=5, » II nj| Ir 
k=n+1 I=n+41 
1Zk 
2r n n+2r 
=> Aid IL HH h (5.5) 
k=1 I=1 l=n+41 
lZn+k 
n 2r n+2r 
= [115 Me H hi 
l=1 k=1 =n+1 
Ps 


Mas por definição, para k = 1, 2,...., 2r— 1 temos 
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Antk fntk41 = Antk4l fntk, 


logo 

n 2r—1 n+2r 
s=[h po ( (+ Dre + 1 Dfre) TJ fi 
Il=1 =1 I=n+1, 

k=9n>1 ltn+k 

Ign+k+1 (5.6) 

n 2r—1 

=[[1 >) 2Re[(1+)fnena] he. 
ERA 


De (5.2) verifica-se que Re [(1 + df .,q] é um polinômio real de grau 1, e de (5.1) 
segue que TE! , é um polinômio real de grau 2n. Portanto, S é um polinômio real de grau 
2n+1+2(1-1) = 2(n+1)-1. Isso prova que P é um polinômio real de grau < 2(n+n-1. 

Resta mostrar que o conjunto de órbitas periódicas de X é igual a C e que cada 
elemento de C é um ciclo limite de X. Para isso introduziremos funções que permitam o uso 
de Teorema 5.3.4. 


Consideremos agora a função definida em Rº — O por 
Hs 


To An+1 An+2r—1 pintar = o 
H == ne O n+1 Do pd À n+2r 


Mostraremos que Hé real e dada por 


H = A(x, ny B(x, yC(x, y), (x, y) E R2 e Q, 


onde n 
As p= ix-y+(y-?-") (5.7) 
1=1 
os) fic-+ (xy) (5.8) 
j=1 


C (x, p=oo(-2) Dr ES =), (5.9) 
5) 


De (5.1) segue imediatamente que o produto dos primeiros n fatores deH 


exatamente 4. 
Para a justificativa de Be C, consideremos 1 <j<re(x, ye R2- (2. Defina Z= 
x— x, + i(y— y). De (5.2), 
Epi 1 AA 2 2 ; j ; + 
Inaja (2,0) fnsoz (2,9) = o PAC) (Oy yr, 
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logo de (5.4),temos 


n+2j—1 n+23 = 
EAR (x,y) feto y) = [(x — 2;)? +ig= y) 2; 2; : 


Assim, tomando A como em (5.7), resulta 


Portanto, tomando B como em (5.8), temos 


H(x,9) = A(e,y)B(e, oTI (2). 


Podemos escrever z com reº, com r, >0. Então, para (x, y) E RP — O x X, temos 


UNE io nº o Nº 

2 roe e QB 
=>) = 5) = =) = (e), 
Z; roeTiê e-iê 
j 


como x < x, temos que Ea logo pela definição de potência complexa 


25) — o, (5.10) 
Zj 


Se x=x, então O = ae 


CRC Eid 


logo (5.10) é válido para x = x, 


De (5.4) e (5.10) segue que podemos escrever 


H (x, y) = A(x, y).B(x, y).C(x, y), 


onde C é dado por (5.9). Assim, a função H é claramente real em R2- 1. 
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Mostraremos agora, que a função V definida por V = | Has f, é fator integrante 


inverso do campo vetorial polinomial X. 
Observemos que da definição de Ve de (5.4) 


V (x 39) = A(x, 9) | (x x) + (y — y)?), (x, ) E RZ (5.11) 


Logo V é uma função polinomial real definida em todo R2. 
Embora Be C sejam positivas, a função A pode assumir valores negativos, logo o 
mesmo ocorre com H . Assim, tomamos R > suficientemente grande de modo que C c 
B(0, R], j=1,2,..., n. Então H >0 no aberto O de Rº definido por (R2 -O)n(R? -BJO, RJ). 
Defina então H = Log F] , que é uma função real O. 
Provemos agora que 
0H Q 90H —P 
= = 5 (5.12) 
dr V' Oy V 


em O,ou seja, que H é uma função hamiltoniana de X/V em O. 
Observemos que, pela Regra de Leibniz para derivação do produto, 


0H OLogH) 10H 1 => ms (AI p )- Du 
É E Ir “TT a. n+2r kJ k = Topo 
a es Hr Pe ge I=Lk mo th 
e 


o. as de Mal 1) E Op s ' Ba 
n+2r = k d 
id Lã k E=1 


Daí segue a primeira igualdade em (5.12). A segunda é análoga e terá sua 


demonstração omitida. 
De (5.12) obtemos 


oV oV oHOV ódHoV 
Cr v(- Oy dx » Ox x). (5.13) 


Por outro lado, podemos reescrever (5.12) da forma 


9H 0H 
=—V Q=V 5 
e Oy Q dx 


Derivando a primeira igualdade com relação a x, a segunda com relação a y e 


somando os resultados obtemos 


oP OP oHOV 0H09V 
+ = + —. (5.14) 
dx Oy dy dr Ox Oy 
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De (5.13) e de (5.14) segue que V é um fator integrante inverso de X em O, ou seja, 
que 


p?V OV (5 3) V. (5.15) 


[oh +03 “Adr a Oy 


em O. Como P Qe Vsão polinômios, a igualdade (5.15) é válida em Rº. 
Pelo Teorema 5.3.4 os ciclos limites de X, se existirem, estão contidos no conjunto 
V=0. De (5.11) resulta 


V(z,y)=06 (2,9) e(Uc]) LI frases pr) (5.16) 


Assim, temos que se X possui ciclos limites, eles devem ser os círculos C, 1 <i<n. 
Resta mostrar que os círculos C, são ciclos limites. 

Mostraremos agora que sobre cada C, não existe um ponto singular de X. Para isso, 
suponhamos que (X,, Y,) seja um ponto singular de Xem C, para algum i= 1,2,..., n, isto é 


P (Xy Yo) = AX, Yo) = HXo Yo) = 0, 


logo, da definição de Pe Q 


n+2r of, 
Plxo,yo) -—( FE tico, no) = (20,0) = O 


l=1 fi 


n+2r 
Q(xo, yo) = o FE tico, no) == (x0,0) = 0. 


l=1 li 


Como f(x, Y)) £ O, para todo |  i, obtemos 


of q Of 
dy º da 


= (. 


Assim, 2(x,- x)=0 e 2(y,- Y) = 0, O que implica que x, = x,e Y, = Y, Portanto, (x, 
Yo) é um ponto do centro do círculo C, o que é uma contradição. Então o círculo C, é uma 
órbita de X. 

Para concluir a demonstração, resta mostrar que cada C, i= 1,2, ..., é um ciclo 
limite. 


Suponhamos que C, não seja ciclo limite. Então existe uma órbita periódica 
v=tx(d, vtd), te R$, 
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diferente de C, 1 sis n, e suficientemente próxima de C tal que em Inty existem os 
mesmos pontos de (p,, ..., p) do que em IntC. Sem perda de generalidade, assuma que 
esses pontos são (p,, ..., pj. 
Note que todos os círculos C,, C,, ..., C, e todos os pontos p,, p,, ..., p, estão na 
curva H=0€ que estas são as únicas órbitas de X nesse nível. Como y é diferente de C,, 
» C, existe h O tal que 


E (xd, vd) = Alx(D, AD) BUD, y(D) exp (=; 204 : Hen 


onde 


g;(t) = arg (ot) 


v(t) — x; 


Notemos que t » A(x(?, y())B(x(), Y(?)) é limitada pois y é compacto. Pela de- 
finição de 0, os ângulos 9 (2, ..., 0 (7) tendem simultaneamente para +º ou -o quando t 
> +ºº, enquanto os ângulos id: -.., 94?) permanecem limitados quando t > +00, como 
podemos ver na Figura o que é uma contradição. Portanto, C, é ciclo limite para o campo 
vetorial X. Como ie (1, ..., ny é arbitrário, o teorema está provado. 

Observe que, na demonstração do teorema fornecemos uma expressão explícita 
para o campo de vetores polinomial satisfazendo o item (b) do Teorema 5.2.2. Claramente, 
o item (a) segue diretamente do item (b). 

Apresentamos a seguir, um exemplo que ilustra o Teorema 5.2.2, o exemplo pode 
ser encontrado em Y. IPYASHENKO [16]. 

Exemplo 5.4.1 Considere C=(C,, C,), onde 

C=2+y/-1=0€C,=2+y/-4=0. 

Neste caso, temos que somente a curva C, é primária. Defina 

fg )=2+/2-1=0, x )=8+/-4=0, 


Realização de uma configuração de ciclos limites 


58 


Figura 5.2 — Ciclo limite C, e órbita periódica v. 


C, dá 


Fonte: Elaborado pelo autor. 
lo y=x+iyeflx y)=x-iy 
e seja a função 
4 
7] A fÃo pA3 pÀ A 
H(x, y) — [ 1 “fo fa? sx, y) = Lg" (x, 1). 
k=1 


comh =h,=1,h,=1+ieA,=1-—i Assim, 


H (x, 3) = (08 + GÊ + yº — 408 + 9º) (er? arcar, 
De onde vemos que a função H é uma função real. Defina 


4 
H =logH =53  Mlog fi. 


k=1 
Notemos que H é também uma função real. 
Definindo o campo de vetores 
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x o fa 


Alhbtdso = Ml foro SiS 


X=(M(fafat SE = dolbitafo ço E 


0 fa 


HOu(fafato GE + dah af SE 4 (fi fat) E + Moto) ES 


E 


dy” 


que no nosso caso torna-se 
—(6x*y +12x2yº —-20X2y +6)º —20yº —2xX —2y'x-4Xyº +10xXº +10yºx+8y —8x) > 
x 


+(6x + 12xºyº — 20Xº + 6y'x— 20y2x + 2Xºy + 2º + 4yºx2 — 10yxº — 10y + 8x+89) À, 
y 


Concluímos que o retrato de fase possui dois ciclos limites (C, e C,) e um equilíbrio 
na origem. Seu esboço se encontra da Figura 4.2. 


Figura 5.3 — Retrato de fase do sistema do Exemplo 5.4.1. 


>< 


Fonte: Elaborado pelo autor. 


Realização de uma configuração de ciclos limites 


60 


REFERÊNCIAS 


1. F DUMORTIER, J. LLIBRE and J.C. ARTÉS; Qualitative theory of planar differen- tial systems. 
Springer, 2006. 


2. F VERHULST; Nonlinear differential equations and dynamical systems. Epsilon Uitgaven, Utrecht, 
Springer-Verlag Berlin Heidelberg, 1985. 


3. H. GIACOMINI, J. LLIBRE, and M. VIANO; On the nonexistence, existence and uniqueness of limit 
cycles. Nonlinearity, 9:501-516, 1996. 


4. J. A. SANDERS, F. VERHULST and J. MURDOCK; Averaging methods in nonli- near dynamical 
systems. Volume 59 of Applied Math. Sciences, Springer, 2000. 


5. J. ÉCALLE; Introduction aux fonctions analysables et preuve constructive de la conjecture de 
dulac. Actualités Math. Hermann, Paris, 1992. 


6. J. K. HALE; Ordinary differential equations. Krieger Publishin Company, 1980. 


7. J. LLIBRE and G. RODRÍGUEZ; Configurations of limit cycles and planar poly- nomial vector 
fields. J. Differ. Equat. 198:374-380, 2004. 


8. J. MUNKRES; Topology. Prentice Hall, Inc. Upper Saddle River, 1975 


9. J. SOTOMAYOR; Lições de equações diferenciais. volume 11 of [Euclid Project ]. Instituto de 
Matemática Pura e Aplicada, Rio de Janeiro, 1979. 


10. L. PERKO; Differential equations and dynamical systems. Springer-Verlag, New York, 1991. 


11. M. G. SOARES; Cálculo em uma variável complexa. Instituto de Matemática Pura e Aplicada, Rio 
de Janeiro, 1999. 


12. M. 1. ALMUKHAMEDOV; On the construction of a differential equation having given curves as 
limit cycles. Izv. Vyss. Ucebn. Zaved Matematika (in Russian), 44:3-6, 1963. 


13. PV. ARAÚJO; Geometria diferencial. Instituto de Matemática Pura e Aplicada, Rio de Janeiro, 1998. 


14. S. SCHECTER and F. SINGER; A class of Vector Fields on Sº That are Topolo- gically Equivalent 
to Polynomial Vector Fields. J. Differ. Equat. 57:406-435, 1985. 


15. X. HUANG, Y. WANG and A. CHENG; Limit cycles in a cubic predator-prey diffe- rential system. 
J. Korean Math. Soc. No. 4, pp. 829-843, 2006. 


16. Y. IPNASHENKO; Finiteness theorem for limit cycles. volume 94 of Trans. of Math. Monogr. Am. 
Math. Soc. 1991. 


17. C. THOMASSEN; The Jordan-Schonflies theorem and the classification of sur- face. The Amer. 
Math. Month., Vol. 99, No. 2 (Feb., 1992), pp. 116-131. 


Referências 


61 


APÊNDICES 


APÊNDICE A - ESCLARECIMENTOS 


Apresentamos aqui a demonstração do Teorema 4.3.3, o qual utiliza uma versão 
do Teorema de Schonflies para grafos. Para isto, precisamos de algumas definições e 
resultados de Topologia do plano. 
Um arco simples no R” é a imagem de uma aplicação injetora contínua do intervalo 
real [0, 1] em Rº. 
Lema A.0.1 Seja O c Rº um aberto conexo, então quaisquer dois pontos distintos de 
O podem ser ligados por um arco simples em O. 


Lema A.0.2 Seja K c R? compacto, Lc R2umaberoe p : K> L é contínua e 
injetiva, então p é um homeomorfismo sobre sua imagem. 


Lema A.0.3 Seja X= AU B, onde X c Rº, Ae Bsão fechados em X. Sef:A5> Ye 
9:B> Ysão contínuas, Yc R2, e f(x) = g(x) para todo xe An B, então a função h: 
X> Y'definida por h(x) = f(x) sexe Ae h(x) = g(X) sex € B, é contínua. 

Os resultados que seguem podem ser encontrados em C. THOMASSEN [17]. Co- 


meçamos com uma versão do Teorema de Schonflies para uma única curva de Jordan. 


TeoremaA.0.4 Sejam, C' curvas de Jordanno planoeh: C> C'um homeomortismo. 
Então existe um homeomorfismo H : Rº > Rº cuja restrição ao conjunto C conincide com h. 

As próximas definições têm por objetivo estender o Teorema de Schonflies para 
grafos. 

Definição A.0.5 Um grafo G é a união de dois conjuntos disjuntos finitos disjuntos 
V(G) e E(G) (cujos elementos são chamados de vértices e arestas, respectivamente), de 
modo que para cada aresta estão associados dois vértices distintos x e y, chamados de 
extremi- dades da aresta. Denotamos uma aresta cujas extremidades são x e y por aresta 
xy. 

Definição A.0.6 Um caminho em G é um grafo com distintos vértices v, ..., v e 
arestas V,V,, ..., V.4V, Sen >2, e quando adicionamos uma aresta vv, para este caminho, 
obtemos um ciclo. 

Se Gé umgrafoe Ac V(G) U E(G), então G — Aé o grafo obtido de G, excluindo 
todos os vértices de 4, e todas as arestas que estão em Aou são incidentes com um vértice 
em 4. Se vé um vértice do grafo G, utilizando um abuso de notação, denotamos por G — Lv) 
o grafo obtido a partir de G, excluindo o vértice ve as arestas que são incidentes com v. 

Definição A.0.7 Dizemos que G é conexo se todos os pares de vértices em G estão 
uni- dos por um arco simples. G é 2-conexo se for conexo e, para cada vértice v, também 


for conexo o grafo G — (v), para qualquer vértice v de G. 
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Definição A.0.8 Um grafo planar é um grafo cujos vértices são pontos do plano, 
as arestas são arcos simples em Rº e além disso, quando duas arestas distintas têm no 
máximo suas extremidades em comum. 

Se G é um grafo planar 2-conexo com pelo menos três vértices, então demonstra- 
se que R2 — Gé a união de IV (G)l + IE(G)l + 2 abertos conexos cada um dos quais tem 
um ciclo de G como fronteira. Cada uma delas é chamada uma face de Ge a única face 
não limitada é chamada de face exterior de G. Aqui os símbolos IV (G)l e IE(G)l denotam a 
cardinalidade dos conjuntos de vértices e arcos de G, respectivamente. A fronteira da face 
exterior é chamada de ciclo exterior. 

Definição A.0.9 Um isomorfismo entre os grafos G e H é uma bijeção entre os 
conjuntos de vértices de Ge H, f: V(G) > V(H), de tal forma que quaisquer dois vértices 
uevde G são adjacentes em G se, e somente se, f(u) e f(v) são adjacentes em H. Nesse 
caso G e H são chamados de grafos isomorfos. 

Definição A.0.10 Um isomorfismo entre os grafos Fe F' é chamado de isomorfismo- 
planar quando, um ciclo é a fronteira de uma face de F se e somente se, a imagem do 
ciclo é a fronteira de uma face de F' e, além disso, o ciclo exterior de F é aplicado no ciclo 
exterior de F'. 

O próximo resultado é uma versão para grafos do Teorema de Schontflies A.0.4. 


TeoremaA.0.11 SejaFeF' grafos planares 2-conexos tal que g é um homeomorfismo 
e um isomorfismo-planar de F em F'. Então g pode ser estendido para um homeomortismo 
sobrejetivo definido em todo o plano. 

Finalmente, passamos agora à prova do Teorema 4.3.8. 

Demonstração: Seja C=(C,, C,, ..., C) a configuração de curvas fechadas dada 
no Teorema 4.3.3. Consideremos o caso n >2, pois o caso n = 1 segue imediatamente do 
Teorema A.0.4. Suponha inicialmente que todas as curvas de C são primárias, tal como 
ilustrado na Figura A.1. 

Para cada j = 1, 2, ..., n escolha um círculo C/de modo que todos os elementos da 
configuração de curvas fechadas dadas por C'=(C,, C,, ..., C ) sejam primárias. Considere 
ainda para cada j= 1, 2, ..., num homeomofismo h,: C, > CG; escolha p, q,€ C,e [op q,€ C, 
de modo que p, + q,e p* Gs 

Utilizando o Lema A.0.1, escolhemos um arco simples g, : [0, 1] > G, c Rº ligando 
p, a p, de modo que a imagem de g, não intercepta as curvas C,, C,, ..., C. Analogamente, 
para cada j = 2, ..., n— 1 escolhemos um arco simples 9;: [0,1] > G, c Rº ligando Ga Ps 
de modo que a imagem de 9; não intercepta as curvas C,, C,, ..., C, tampouco as imagens 
das curvas g,com 1 <i<j. Consideramos ainda 9, : [0, 1] > G, c R? como um arco simples 
ligando q, a q, de modo que a imagem de g, não intercepta as curvas C,, C,, ..., C, 
tampouco as imagens das curvas g, 1 si<n— 1. Analogamente definimos 9; para j=1, 
Bino 1,80 
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Do Lema A.0.2 decorre que cada 9; é um homeomorfismo sobre sua imagem G, 

Análogo para g,. 

Denotamos por Fo grafo 2-planar cujos vértices são os pontos p,, P,. ... PQ» Gy 
«- q, e cujas arestas são os elementos de C bem como os arcos de cada G; ligando p, a q, 
j=1,2,...., n. Analogamente definimos o grafo F. 

Definimos então a aplicação h: F > F' do seguinte modo. Se x € C, para algum j= 1, 
2,..., n, então h(x) = ht), esexe G, então h(X) = J; (a; (a). Pelo Lema A.0.3 segue que 
h é um homeomorfismo, e da construção segue ainda que h é um isomorfismo planar. Do 
Teorema A.0.9 resulta que as configurações de curvas fechadas Ce C' são equivalentes. 

Suponha agora que exista um elemento de C que não é primária, digamos C,, C,, 
«oC CIntC eC,. é curva primária. Nesse caso ligamos as curvas C,, C,, .... C entre 
sital como se C,, C,, ..., C fossem primárias, e depois ligamos um ponto de C, a um ponto 
de C,., e um ponto de C a outro ponto distinto de C 


met? 


como mostra a Figura A.2. 


Figura A.1 — Configurações topologicamente equivalentes. 


c 


Fonte: Elaborado pelo autor. 
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Figura A.2 — Configurações topologicamente equivalentes. 


Fonte: Elaborado pelo autor. 
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